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ABSTRACT 

These lectures, given at the 1997 TASI Summer School, describe the 
prospects for discovering supersymmetry (SUSY) and for studying its prop- 
erties at the Large Hadron Collider (LHC) at CERN. If SUSY exists at a 
mass scale less than 1-2 TeV, then it should be easy to observe characteristic 
deviations from the Standard Model at the LHC. It is more difficult to deter- 
mine SUSY masses because in most models there are two missing particles Xi 
in every event. However, it is possible to use various kinematic distributions 
to make precision measurements of combinations of SUSY masses and other 
quantities related to SUSY physics. In favorable cases such measurements at 
the LHC can determine the parameters of the underlying SUSY model with 
good accuracy. 
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Supersymmetry Signatures at the CERN LHC 



Frank E. Paige 
Physics Department 
Brookhaven National Laboratory 
Upton, NY 11973 USA 

These lectures, given at the 1997 TASI Summer School, describe the prospects for 
discovering supersymmetry (SUSY) and for studying its properties at the Large 
Hadron Collider (LHC) at CERN. If SUSY exists at a mass scale less than 1-2 TeV, 
then it should be easy to observe characteristic deviations from the Standard Model 
at the LHC. It is more difficult to determine SUSY masses because in most models 
there are two missing particles X\ i n every event. However, it is possible to use 
various kinematic distributions to make precision measurements of combinations 
of SUSY masses and other quantities related to SUSY physics. In favorable cases 
such measurements at the LHC can determine the parameters of the underlying 
SUSY model with good accuracy. 

1 Introduction 

The theoretical attractiveness of having supersymmetry (SUSY) at the elec- 
troweak scale has been discussed by many authors .Era But while SUSY is per- 
haps the most promising ideas for physics beyond the Standard Model, we 
will not know if it is a correct idea until SUSY particles are discovered ex- 
perimentally. LEP might still discover a SUSY particle, but it has already 
run at y/s — 183 GeV, and its reach will be limited by its maximum energy, 
probably \fs = 193 GeV. It is less unlikely that LEP might discover a light 
Higgs boson: the current boundB of ~ 77 GeV is expected to be improved to 
& 95 GeV,Q whereas the upper limit on the mass is ~ 130 GeV in the Minimal 
Supersymmetric Standard ModeB and ~ 150 GeV more generally.Q Finding a 
light Higgs would not prove the existence of SUSY, but it would certainly be 
a strong hint. 

The TevatronJias a better chance of finding SUSY particles, particularly 
from the processor] 

pp -» xlxt -> tr^flr ■ 

This can be sensitive to Mfjfy) ~ M(\f) 200 GeV for some choices of the 
other parameters, e.g., small tanJ3, given an integrated luminosity of 2fb _ in 
Run 2 and more in future runs.EI But, like LEP, the Tevatron cannot exclude 
SUSY at the weak scale. 

The decisive test of weak scale SUSY, therefore, must await the Large 
Hadron Collider (LHC) at CERN. The LHC can detect gluinos and squarks in 
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the MSSM up to ~ 2 TeV with only 10% of its design integrated luminosity per 
year. Discovering gluinos and squarks in the expected mass range, ^ 1 TeV, 
seems straightforward, since the rates are large and the signals are easy to 
separate from Standard Model backgrounds. Other SUSY particles can be 
found from the decays of gluinos and squarks. The difficult problem is not 
discovering SUSY if it exists but verifying that the new physics is indeed SUSY, 
separating the various SUSY signals, and interpreting them in terms of the 
parameters of an underlying SUSY model. This is more difficult for the LHC 
than for an e + e~ machine of sufficient energy, but some progress has been 
made recently. 

The first few sections of these lectures are mainly review. Section || re- 
views the SUSY production cross sections and some basic facts about QCD 
perturbation theory. Section || discusses event generators, which are used to 
translate production cross sections into experimental signals. Section ^ sum- 
marizes the capabilities of ATLAS and CMS, the two main LHC detectors, to 
detect these signals. 

The next sections concentrate on SUSY measurements at the LHC in the 
context of the minimal supergravity (SUGRA) model,£3 although the general 
results should apply to other models, at least those in which the lightest SUSY 
particle escapes the detector. Section || shows the reach in SUGRA parameter 
space for various signals and describes how to make a first estimate of the SUSY 
mass scale. Section ^ describes examples of a recently developed approach to 
extracting information about SUSY masses and other parameters from LHC 
data for five specific SUGRA points. Section ^ shows what the resulting errors 
on the SUGRA parameters would be at these points. Section || discusses 
preliminary results at a SUGRA point with large tan (3 that has very different 
properties. 

Section || discusses the LHC discovery potential for Higgs bosons both 
in the Standard Model and in the MSSM. Finally, Section [h] draws some 
conclusions. 

2 SUSY Cross Sections 

The LHC is a pp collider to be built in the existing LEP tunnel at CERN 
(the European Laboratory for Particle Physics, located near Geneva, Switzer- 
land) with a center-of-mass energy y/s = 14 TeV and a luminosity C = |1Q 33 ~ 
10 34 cm~ 2 s _1 . It will have two major experiments, ATLAaB and CMS,E3 for 
studying high-px physics like SUSY. Two smaller experiments — LHC-B for 
B physics and ALICE and for heavy ion physics — have also been proposed 
but will not be discussed further here. Construction of both the accelerator 
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and the experiments is expected to be completed in 2005. El 

Fine tuning argumentsEj suggest that the SUSY masses should be below 
about 1 TeV if SUSY is relevant to electroweak physics. Then SUSY production 
at the LHC is dominated by the production of gluinos and squarks. The 
elementary g and q cross sections only depend on the color representations 
and spins of these particles — which of course are fixed by supersymmetry 
- and on their masses. Thus they are less model dependent than the cross 
sections for gaugino production, which also depend on couplings determined 
by the mixing matrices. 

Perturbative QCD tells us that inclusive production cross sections can be 
computed as a power series in the strong coupling coupling a s (Q) evaluated 
at a scale Q of order the masses involved. For examplc_.the lowest order 
contribution to the elementary process gg — > gg is given byt-51 



where s, t, u are usual parton process invariants and M = Mg is the only 
relevant mass in this case. The lowest order cross section for gg — > qq depends 
on both Mq and M g . The cross sections for gaugino pair production and 
associated production depend on the masses and also on couplings which are 
determined by the chargino and neutralino mixing matrices. 

The elementary cross sections are then related to pp cross sections by the 
QCD-improved parton model, which is based on the impulse approximation for 
processes with large Q 2 . In the parton model the pp cross section is given by a 
convolution of the elementary parton-parton cross sections and the appropriate 
parton distributions, i.e., the probabilities of finding quarks or gluons with 
given momentum fractions Xi in the incoming protons: 



Here Q 2 is a measure of the scale, e.g. Q 2 = p^+M 2 , and Xi are the momentum 
fractions of the incoming partons. By elementary kinematics they satisfy 




(j - M 2 )(u — M 2 ) — 2M 2 (t + M 2 

(t-M 2 ) 2 
M 2 {s-AM 2 ) \ 



<-> u) + 



{t-M 2 )(u-M 2 ) J ' 




s = 
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Figure 1: g. u, and— m parton distributions xf(x,Q 2 ) in the proton vs. x for Q = 100 GeV 
using the CTEQ3LE3 parameterization. 




P 



Figure 2: Soft and collinear singularities both come from gluons attached to external lines, 
only one of which is shown. Soft singularities come from p 2 = m 2 and k — » 0; collinear ones 
come from p 2 = k 2 = and 9 V ^ — > 0. The gluon can be either real or virtual. 



where y is the rapidity of the center of mass of the produced system. 

Representative parton distributions for Q — 100 GeV are shown in Fig- 



so that g(x) is large and gg processes dominate the production of squarks and 
gluinos for most masses of interest. 

The QCD-improved parton model is intuitively plausible, but we must 
require that it is consistent with higher order perturbation theory. Consider 
adding one gluon emission or loop. Then after renormalization of ultraviolet 
divergences in the usual way, individual graphs are found to give a series not 
in a s but in a s \n 2 Q 2 . These logarithms must cancel if perturbation theory 
is to be usable. There is only one large Q 2 in the problem, so the logarithms 
must reflect singularities as the external particles are put on mass shell. There 
are actually two distinct types of singularities, both well known from QED: 



ure 



|. Note that 




Soft or Infrared Singularities: These arise from a gluon with k — > at- 
tached to an external line, Figure ||. Then as the external line goes on shell, so 
does the internal propagator. Soft singularities arise from the divergent multi- 
plicity of soft gluons; the probability of radiating no gluons from an accelerated 
color charge is zero. The total cross section must be finite — something must 
happen. Hence the singularities must cancel between processes with different 
numbers of gluons, i.e., between real and virtual graphs This can be proven 
to all orders in perturbation theory for QED.E3E3 The situation in QCD is 
more complicated because the gluons themselves radiate, but the cancelation 
certainly works in all cases that have been tried. 

Collinear or Mass Singularities: These arise from the emission of a hard 
gluon parallel to a massless parton in the initial state, i.e., p 2 = k 2 = 
and k\\p in Figure || . (There also can be collinear singularities in the final 
state for differential cross sections.) They do not cancel for initial states like 
hadrons with limited transverse momenta, but they are universal because they 
come from on-shell poles. tj Hence they cancel if you calculate one physical 
process, e.g., gg production, in terms of another, e.g., deep inelastic scattering. 
Equivalently, they can be absorbed into universal parton distributions defined 
by deep inelastic scattering. It is straightforward to verify this at one loop. The 
general proof is complex because soft and collinear singularities get tangled,! 2 !! 
but the result is believed to be true. 

The collinear singularities in the parton distributions lead to a series in 
(a s (Q 2 ) InQ 2 ) = 0(1). These leading logarithms can be summed to all orders 
in perturbation theory: Altarelli and ParisiEil matched th^ operator product 
expansion to the 1-loop graphs, while Gribov and LipatovcJ and DokshitserEa 
studied the origin of the logarithms in perturbation theory directly. The result 
is that the parton distributions satisfy the DGLAP equations, 



dq l (x 1 Q 2 ) 
dQ 2 



dg(x,Q 2 ) 
dQ 2 



2ttQ 2 



a s {Q 2 ) 
2ttQ 2 



J^[ qi{x \Q 2 )P qq (^)+ g ( x >,Q 2 ) Pqg (^) , 



where the DGLAP functions 




PuM = ^(x 2 + (l-x) 2 ) , 
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P gg (x) = 2c A [j—^ + + x 0- ~ x )j + CggSi 1 - x ) > 

reflect the various QCD couplings. The + subscript in these functions indicates 
that the singularity is to be treated as a distribution, 



(l-x)+ J (1-z) " 

The 1/(1 — x) singularities come from the radiation of soft gluons, while the 
(5(1 — x) terms come from virtual graphs. A little thought will show that using 
such a distribution implements the real-virtual cancellation just like the usual 
perturbative calculation. It is not actually necessary to calculate the virtual 
graphs: the coefficients of the delta functions can be determined by momentum 
conservation^ 

The DGLAP equations correspond to a simple picture of the Q 2 evolu- 
tion of parton distributions. As Q 2 increases, more gluons are radiated, so 
the distributions soften at large x and increase at small x. This radiation is 
responsible for the rise at small x in Figuxciy. 

The lowest-order QCD cross sectionscSQ for SUSY particle production in 
the MSSM are shown in Figure These cross sections are similar but not 
identical to those in the SUGRA model; the gaugino masses are scaled from g 
like 

Mi _ M 2 _ M g 

Oil «2 Oi s 

instead of being calculated from the renormalization group equations. How- 
ever, the gluino and squark cross sections are model independent, and these 
dominate except for very high masses and heavy squarks. 

In addition to their effects on parton distributions, higher order QCD cor- 
rections also give finite 0{a s ) corrections to pp cross sections. These correc- 
tions typically give significant corrections to the overall normalization. That 
is, they give 

a k, Kg lo 

with a "if -factor" 

7T 

that is typically substantially larger than one but less than two for the natural 
lowest-order scale choice because C ~ ir 2 rather than C ~ 1. While the 
effect on the normalization is significant, the effect on the shape of inclusive 
distributions is typically small. This is illustrated in Figure ^, which shows the 
Pt distributions for the leading-order (LO) and next-to-leading-order (NLO) 
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Figure 3: Total SUSY MSSM production cross sections'] at the LHC using M; = otiMg/a s 
to determine the gaugino masses. 



calculations for SUSY production at the LHC. The two shapes are clearly 
almost identical. 

The similarity in the shapes of the LO and NLO distributions seems to 
be a general feature of perturbative QCD calculations for a wide variety of 
pp proposes. It is partly understood, at least for the simple case of Drell- 
Yan.t3E3 The NLO Drell-Yan calculation contains a ln 2 (— Q 2 ) factor from the 
overlapping soft and collinear singularities. The logarithm is canceled by a 
ln 2 (Q 2 ) factor for deep inelastic scattering, but there is a tt 2 from the difference. 
This tt 2 multiplies the natural scale of (4a s /3ir)ao, where <jq is the lowest order 
cross section, so it produces an overall normalization factor but no change in 
shape. It is not clear, however, that this 7r 2 factor is the dominant effect even 
in the simple case of Drell-Yan. 

The cross sections shown in Figure |^ correspond to large rates. "Low" lu- 
minosity at the LHC is 10 33 cm~ 2 s (about 10 2 times that currently achieved 
at the Tevatron) while the design luminosity is 10 34 cm _2 s _1 . A "Snowmass 
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Figure 4: Normalized LO and NLO pr distribution for squarks and gluinos at the LHC 



Table 1: Events/year for various SUSY masses at LHC at low and high luminosity using the 
cross sections shown in Figure pi 
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year" is defined to be 10 7 s, not 3.1536 x 10 7 s; it represents the typical ef- 
fective running time per year for an accelerator. For example, in the recently 
completed Tevatron Run I, an integrated luminosity of about 100 pb -1 was ob- 
tained in about two years with peak luminosity of about 10 31 cm s . Given 
a luminosity of 10 4 — 10 5 pb _1 per year, a large number of SUSY events will 
be produced if the masses are below ~ 1 TeV, as can be seen from Table [l]. 

The gluino and squark px distributions must obviously peak at px ~ 
M, since smaller px's are suppressed by phase space, while larger ones are 
suppressed by the requirement of additional energy for the initial partons. This 
peaking can be seen in Figure [l| both for the LO and for the NLO calculations. 
This means that the gluino or squark decay rest frame and the lab frame are 
similar, and the decay products in the lab frame are spread out over phase 
space. Gluinos or squarks decay via several steps to the lightest SUSY particle 
(LSP), taken to be the Xi- Thus a typical event might be: 

9 + 9 -> 9 + 9 



8 



9 -> Xi<7 • 

This example event contains five new particles(l). The Xi has electroweak cou- 
plings and must produce virtual sleptons or squarks when interacting with mat- 
ter, so interacts weakly and escapes from the detector, giving fx- Signatures 
for such events involve multiple jets, large fx from the two weakly-interacting 
Xi's, possible leptons from decays of gauginos or sleptons in the cascade decay 
process, and b jets that are either democratically produced or enhanced by 
smaller third-generation squark masses and large Yukawa couplings. 

All of these features provide handles that can be used to separate the SUSY 
signals from Standard Model backgrounds. As will be discussed in Section || 
below, the LHC detectors can measure the Ex of all jets and leptons, so they 
can determine fx (with possible problems from cracks and resolution tails). 
They can identify and measure electrons and muons well; hadronically decaying 
r's can also be identified, although with more background. Finally, b jets can 
be tagged using vertex detectors to observe displaced vertices. 

The single most important experimental signature for SUSY at the LHC 
is probably missing transverse energy fx- The Standard Model backgrounds 
for fx come from all the possible ways to produce high-p^ ^'s: 

• W — > £v, tv 

• Z — ► VV, TT 

• QQ~* £uX, Q = c,b,t. 

• 9 QQ, Q -> tvX 

The electroweak cross sections are suppressed by powers of a, while the heavy 
quark ones are suppressed by color and spin factors. Both are suppressed 
by leptonic branching ratios and by the fact that the missing v typically has 
smaller px than its parent. There are also backgrounds to fx from cracks and 
resolution tails in the detector. These are difficult to estimate in any simple 
way, but detailed calculations indicate that real neutrinos dominate for events 
with multiple jets and/or leptons and fx i~ 100 GeV. 

3 Event Simulation 

SUSY signatures typically involve multiple jets and/or leptons plus missing 
transverse energy fx arising from complex cascade decays. The Standard 
Model backgrounds for these signatures arise only from high orders of per- 
turbation theory, and the dominant contributions come from those regions of 
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phase space that give soft or collinear singularities. Thus it is appropriate to 
use Monte Carlo event generators that simulate complete events and include 
the most important effects of QCD radiation to calculate both the signals and 
the backgrounds. Such event generators also provide complete events that can 
be used to estimate detector-induced backgrounds. ._. 

Three. general-purpose_event generators are in common use: HERWIGcS, 
ISAJETE3, and PYTHIA.E3 Of these, HERWIG provides the best theoretical 
treatment of QCD but presently contains no SUSY processes. ISAJET has the 
most detailed treatment of SUSY. PYTHIA provides a well-tuned description 
of QCD and hadronization plus many SUSY processes. While a precise descrip- 
tion of QCD and hadronization may eventually be important, it it probably 
not essential for the exploratory studies now being done. 

Any Monte Carlo event generator requires a random number generator. 
If the numbers were truly random, then each run would produce different 
results, and debugging would be nearly impossible. Hence we require a pseudo- 
random number generator that produces the same sequence of random numbers 
given the pune initial seed(s). A commonly used algorithm is the congruential 
generator .till It is surprisingly simple: 

so > 0, 
s n +i = (as n + c) mod to . 

For careful choices of a and c and a large value of to, e.g. a = 7 5 , c = 0, to = 
2 31 - 1 or a = 6909031764870, c = 0, to = 2 48 , this algorithm yields a uniformly 
distributed sequence of integers satisfying most tests of randomness. (It is 
obvious that there are also extremely bad choices.) Then x n — s n /m provides 
a uniformly distributed real number in (0, 1). The congruential generator does 
have the limitation that ^-tuples of x n lie on at most rn}l k planes in fc-space. 
The random number algorithm does not at this time seem to be the critical 
issue for LHC SUSY physics, so we will not discuss it further here. 

3.1 Hard Processes and Perturbative QCD 

Given a random number generator, one can construct an algorithm for gener- 
ating events according to any desired parton processes. This section outlines 
the basic steps, but since this Summer School is devoted to SUSY rather than 
to QC TT it opts for simplicity rather than attempting to review the state of 
the artll 

Step 1: Generate the parton process, i.e., pick the kinematic variables 
according to the appropriate lowest order parton cross section. Using a lowest 
order cross section is important: higher order cross sections would involve 
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singular distributions. These can be treated in principle using a cutoff, but the 
resulting distributions have very large positive/negative weight fluctuations. 

This step is simple in principle, although it may be complicated in practice. 
From the definitions of cumulative probability distributions and random vari- 
ables, if x is distributed according to a (normalized) probability distribution 
f(x), < x < 1, then 



where C is a uniformly distributed random variable that can be generated 
using a congruential or other generator. In a few cases, e.g. f(x) = (n + l)x n , 
this relation can be solved analytically. Generally, an analytic solution is not 
possible, but one can find a bound f(x) < / max and use a rejection algorithm, 
accepting x = £1 if f(x) > £2 /max- It is possible to combine these methods 
and/or change variables to improve the efficiency. 

To generate a parton process like gg — > gg, one should choose a set of 
variables like PtiV\tV2i < \> to isolate the rapid cross section variation in a just a 
few variables, e.g., pt- For processes like gluon and light quark jets that do not 
involve a mass scale, some form of importance sampling or weighting should 
be used to generate px efficiently. For process like Z or H with a narrow s- 
channel resonance, one should instead choose a set of variables such as s, i, y, <f). 
Generating the hard scattering is basically straightforward, although it can 
involve rather complex computer code. 

Step 2: Add QCD radiation. Perturbation theory describes inclusive cross 
sections, but it does not necessarily give a good description of event structure. 
QCD is approximately scale invariant, with a dimcnsionless coupling a s whose 
running with Q 2 is only logarithmic. Hence, radiation at all Q 2 scales is 
important; it produces fixed angle jets with <C (M 2 ct ) ~ a s (p 2 ^)p 2 ^ -C p\. 

The most important QCD radiation is the emission of gluons collinear 
with initial or final partons. But we know from the factorization theorems of 
perturbative QCD that collinear singularities factorize, not for amplitudes but 
for cross sections. That is, 



where Pjj is the appropriate Altarclli-Parisi (DGLAP) function and z 
The starting point for QCD-based event generators like HERWIG, ISAJET, 
and PYTHIA is to use this collinear approximation for gluon radiation to- 
gether with exact, non-collinear kinematics. This is known as the branching 




2irp 2 



Pij(z), P 2 ^0, 
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approximation. The branching approximation correctly describes the dom- 
inant, leading-log QCD effects, but because it uses exact kinematics, it also 
gives a fairly good approximation for real higher-order QCD effects such as 
multiple jet production. 

The discussion that follows is limited only to the simplest case of leading- 
log gluon radiation from an outgoing quark line. This does not represent the 
state of the art, but it is simple to present and illustrates the main point, 
namely that a classical branching process can provided a reasonable approxi- 
mation to QCD to all orders in perturbation theory. The cross section for the 
radiation of n gluons from an outgoing quark is given in the collinear limit by 



1 



da 



a dp\dz\ . . . dp„dz 7 , 



as(pl) 
2irpj 



2-kvI 



-P(Zn) 



where z is the momentum fraction. While z is well defined in the collinear 
limit, its non-collinear extension is model dependent. A good choice is 



E, 



+1 



\p 



1+1 1 



qq is 



Ei + \m\ 

With this choice, radiation from the quark and antiquark in e + e~ 
treated symmetrically. 

Introduce a minimum mass t c to regulate both the soft and the collinear 
singularities. For a branching to — > t\ + t 2 in the initial to rest frame, the final 
energy and momenta are given by simple two-body kinematics: 



E* = 



to + t\ — t 2 



Pi 



[(to - h - t 2 f - Ahh} 1 / 2 
2jto~ 



The minimum energy and momentum of decay product 1 depend on the rela- 
tion between the mass and velocity of the decaying system, but the maximum 
energy and momentum always correspond to cos 8—1. Boosting the center of 
mass momenta back to the lab frame, we find 



= 1 



t +h-t 2 , [(t -h -t 2 f -At x t 2 Y' 2 



2to 



2tn 



Setting t\=t 2 = t c , the cutoff mass, gives 



1 - 



4ic 
to 



This cutoff will be used in what follows. 
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The next step is to sum up the radiation softer than the cutoff into a 
Sudakov form factor, i.e., into the probability of evolving from an initial mass 
t to a final mass t' emitting no resolved radiation greater than the cutoff. The 
calculation of this Sudakov form factor is particularly simple if one takes the 
cutoff to be not a fixed mass t c but a fixed smallest z c defined by the initial 
mass and t c cutoff. Kinematics requires 



By definition of a z cutoff, all the Zi are in the interval 

K = {z < z c } U {z > 1 - z c } . 

Thus, using the factorized form for gluon emission, the Sudakov form factor S 
is given by 



The form of this expression is simple because using a z cutoff decouples the z 
integrals from the t ones. 

Virtual graphs contribute 6(1 — z) terms to the DGLAP functions. Mo- 
mentum conservation requires that the ^-weighted integral of the gluon and 
quark distributions give the total momentum fraction, unity: 



This fixes the constants c qq and c gg in the DGLAP functions to satisfy 



where the 1/(1 — z) singularities in the functions P(z) are regulated using "+" 
distributions: 



The physical meaning of these + distributions is that the soft singularities 
cancel between the real and virtual graphs. Because the complete integral 



t > h > . . . > t' . 
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over P vanishes, the z integrals in the definition of S over the unresolved soft 
radiation in the region 1Z can be related to a well-defined integral over non-soft 
values of z in the complement of 1Z: 

/ dzPi(z) = - dz Pi(z) = --fi(z c ) . 

Jtz Jz c 

The integrals in the definition of S over the U are elementary, and the nesting 
gives a factor of 1/n!. The final result has the simple form 



&7«(*c) 

Si(t,t') 



a s (t) 



a s {t>) 

Knowing the cumulative distribution S is enough to generate the next 
mass and hence the complete shower. Given S, the probability E(t) that the 
first resolved radiation occurs at t is given by 

f dtZ{t) = 1-Si(t ,t c ). 

The mass-squared t of the quark at the next branching is therefore generated 
by solving the equation 

1 — Si (t, t c 



1 - Si(t ,t c ) 



C, 



where £ is a uniformly distributed random number in (0, 1). Given the form 
for S, this equation can be easily solved in terms of elementary functions. This 
solution can then be used to build the complete shower: 

1. Generate t. 

2. Generate z c < z < 1 — z c according to Pi. 

3. Check z against t c limits. 

4. Generate t\, ti starting at z, (1 — z)t. 

5. Solve 2-body kinematics. 

6. Iterate. 

This simple Monte Carlo algorithm generates all the leading-log effects of QCD 
radiation and provides a reasonable approximation to non-leading effects such 
as multi-jet production. 
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The algorithm just described can be extended to include g — > gg, qq and 
initial state radiation. It can also be refined in various ways, in particular to 
include in an approximate way coherent interference effects that reduce the 
radiation of soft gluons. For a recent review see Webber .l3 

Step 3: Generate decays of the SUSY particles. This involves selecting 
the decay mode using branching ratios calculated from the SUSY masses and 
couplings and then then generating the decay using 2-body or 3-body phase 
space. Parton showers from any outgoing quarks or gluons are added using 
the branching approximation with an initial scale equal to the parent particle 
mass. 

The use of phase space for a 3-body decay, say q — > x\qq, is not a good 
approximation if the squark q is just slightly heavier than the gluino. The 
squark poles should cause the matrix element to peak where the gaugino and 
one quark are hard, while the other quark is soft. All distributions should be 
continuous as the squark mass varies from below to above the gluino mass, 
but the phase space approximation is not. This problem does not affect any 
of the points to be studied here. The approximations also lose information on 
spin correlations between particles, but this is probably not very important, 
at least for LHC studies. 

Step 4: Hadronize the partons, i.e., generate a jet of hadrons for each 
parton. This involves soft, nonperturbative QCD physics, so one must rely on 
models tuned to experimental fragmentation data. The physical picture is that 
the hard scattering creates separated color charges connected by "strings" of 
gluons. As the strings stretch, they break, pulling a qq pair out of the vacuum. 
Since QCD is strong only at low Q 2 , this breaking is a soft process which 
only occurs with low pr and small Ay. Thus it produces a jet of hadrons 
with limited pr relative to the parton direction. Hence, the event structure is 
mainly controlled by perturbation theory, not by the soft physics model. 

The hard part of a pp interaction never accounts for more than a small 
fraction of the total energy. The rest of the energy is carried off by spectators, 
which produce low-px hadrons more or less uniformly in rapidity. Typically 
the leading particle in the beam jet carries half of the beam energy and just a 
few hundred MeV of transverse momentum. As a result, much of the energy 
escapes down the beam pipe: one can measure missing transverse energy 
but not missing total energy. The models for beam jets are essentially pa- 
rameterizations of experimental data. Fortunately, the beam jets are not very 
important for the experimental signatures, although they do have some effect 
on the isolation requirements for leptons and photons. 
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3.2 Soft or Ins Physics 

The production of SUSY particles and of all the significant Standard Model 
backgrounds for them are hard processes that can be calculated in QCD per- 
turbation theory. But most of the events at the LHC are soft interactions, for 
which there is no good theoretical description. These events are not important 
as backgrounds for new physics, but they dominate the rate and so are impor- 
tant for the design of the detectors. Fortunately, soft physics varies slowly with 
Ins, so reliable extrapolations from lower energy can be made even without a 
good theoretical model. 

The geometrical size of the proton is set by 1/M W , implying a natural 
geometrical cross section 

(he \ 2 
I w 62 mb 
Mj 

since he = 197MeV-f. This should be compared with SUSY cross sections 
ranging from 0.01 pb to 100 pb depending on the masses. The total cross section 
is of course given by the imaginary part of the forward scattering amplitude, 

o-t = -3/(a,0). 
s 

In the high energy limit with fixed momentum transfer (s — > oo, t fixed) the 
exchange of an elementary particle of spin J in the t channel leads to 

f( Sl t)~Pj(z t )~s J . 

The Froissart bound, based only on analyticity properties that can be proven 
in field theory and on unitarity of the S'-matrix, limits the growth of the cross 
section to 

o~t < const In 2 s 

so we must have J < 1. 

The experimental data on the pp and pp total cross sections, Figure ^, show 
rapid variation at low energy and large differences between the pp and pp. This 
difference is not surprising, since pp has a large cross section to annihilate into 
mesons, while pp cannot do this. At high energy, however, the pp and pp cross 
sections become equal and vary slowly with In s, in a way quite consistent with 
the Froissart bound. 

While the total cross section is certainly nonperturbative, it is useful to 
consider how multiparticle production works in perturbation theory. Because 
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Figure 5: Compilation of pp and pp cross sections. 




Figure 6: The discontinuity of the box graph gives a constant cross section for s — » oo. 



the gluon has spin one, the box graph, Figure [| gives a constant cross sec- 
tion for large s. Of course the cross section is infinite for a massless gluon 
exchanged between unconfined quarks, so an infrared cutoff is needed. Hence 
this whole discussion is only qualitative. However, it is possible to do a rigor- 
ous calculation of the behavior of the high-p^ jet cross section for s ^> jij, in 
QCD perturbation theoryEj 

The generalized ladder graph in Figure ^ similarly remains constant as 
the generalized rungs, the shaded blobs in the figure, are moved in rapidity. 
Hence, this graph gives a factor of Ins from each longitudinal phase space 
integral. It can be shown that graphs like this give the leading powers of In s 
for each order in perturbation theory. The n longitudinal integrals are ordered 



17 



Figure 7: The discontinuity of the (generalized) ladder graph gives a factor of In s from each 
longitudinal phase space integral. 



in rapidity and so give 



In s py 



JO 



1 



dyi... dy n = — ln n s . 



ii 



The discontinuity across the generalized rungs produces particles in the final 
state, and the distribution of these particles, like that of the rungs, is flat 
in rapidity. Furthermore, after many steps in rapidity, any quantum number 
exchange will be washed out, so the distribution in the central region is univer- 
sal. All of these properties of this toy model are in qualitative agreement with 
experiment. However, the growth of Figure ^ with Ins violates the Froissart 
bound, so more complicated graphs must become important. 
Data on the charged multiplicity dN/dr), where 

E+ PL 

77 = — In tan — w y = in tor px 2> 777 , 

2 E - p L 

are shown in Figure ^| for pp and pp interactions. The rise is consistent with 
either a power of In s or a small power of s. A smooth extrapolation to 14 TeV 
gives 

dN ch _ 
dr/ 

The mean px also grows slowly with y/s, rising from about 0.35 GeV at low 
energy to about 0.5 GeV at 1.8 TeV. Most of this rise probably comes from 
QCD jets with transverse momenta of a few GeV. Extrapolation suggests 
(p T ) ~ 0.65 GeV at LHC. 



18 




Soft interactions do not produce backgrounds for SUSY or other interesting 
processes, but they do give a very high interaction rate. The LHC bunch 
crossing rate is 40 MHz, and there are ~ 2 interactions/bunch at 10 33 cm _2 s _1 
and ~ 20 interactions/bunch at 10 34 cm~ 2 s _1 . 



4 LHC Detectors 

Two LHC large detectors, ATLAS0 and CMsJll are just beginning construc- 
tion. Both ATLAS and CMS are designed to identify and measure all the 
Standard Model quanta — 7, e, /i, r, g and q jets, and b jets — over \rj\ 2.5 
(6 > 10°) and to measure energy flow over \rj\ ^ 5 (9 > 1°) to determine the 
missing transverse energy fx- This broad coverage is essential for complex 
signatures like SUSY. The cost of each detector by CERN accounting rules is 
CHF 475M.|] Each collaboration currently includes about 1500 physicists and 
senior engineers. Thus, both ATLAS and CMS are more like laboratories than 
traditional experiments. 

The cost and scale of the detectors is driven both by the need to mea- 
sure the high energies and large range of possible interesting processes such 
as SUSY production and by the need to cope with the high rates caused by 
soft physics at the LHC. While the detectors are quite different in detail, they 
both have the same basic elements. These are from inside out: a silicon vertex 
detector intended primarily for tagging b jets, an inner tracker to measure the 

"Both the U.S. and CERN do not include physicists salaries. CERN also does not include 
most EDIA (Engineering, Design, Inspection, and Administration) and labor costs. This 
means that 1 CHF at CERN is roughly equivalent to $2 in the U.S. 
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Figure 9: Cutaway view of the ATLAS detector 11 as of 1997, showing from inside out the 
central detector, 2 Tesla solenoidal magnet coil, liquid argon electromagnetic calorimeter, 
hadron calorimeter and muon system including the lumped toroidal magnet coils and three 
layers of muon chambers. 
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Calorimeters 

Figure 10: Cutaway view of the ATLAS calorimeter and central tracker, the central part of 
Figure g. 



momenta of charged particles in a magnetic field, an electromagnetic calorime- 
ter to measure the energies and directions of photons and electrons, a hadron 
calorimeter to measure jets, a forward calorimeter mainly to measure $t but 
also to tag forward jets, and a muon system to identify and to measure muons. 



These parts can be seen in Figures |9| - |12j. Section 4A below gives a brief 
description of each of these elements and their most important performance 
parameters. This is followed by Section which explains how the parts are 
used in combination to detect physics signatures. The discussion given here is 
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CMS 

A Compact Solenoidal Detector for LHC 




Figure 11: Cutaway view of the CMS detector as of 1997, showing from inside out the 
central tracker, crystal electromagnetic calorimeter, hadron calorimeter, 4 Tesla solenoidal 
magnet coil, and muon system. 
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Figure 12: Cross section of one quarter of the CMS detector, showing the same elements 
as Figure |l l|. 



necessarily superficial, but hopefully it will prove useful to theorists interested, 
in LHC physics. For details the reader should see the ATLASEl and CMSE3 
Technical Proposals and the Technical Design Reports for the detector sub- 
systems. The Particle Data GroupE-3 provides useful general information on 
particle detectors. 



4-1 Detector Elements 

Silicon vertex detector: Both ATLAS and CMS have silicon microstrip detec- 
tors covering \r/\ < 2.5. These are primarily intended to tag b jets by detecting 
the displaced vertices of B hadrons, but they also contribute to the momentum 
measurement. The detector elements are made out of silicon wafers similar to 
those used to make computer chips. A charged particle passing through the 
silicon causes ionization, and the resulting electrons are collected on strips 
about 80 /im wide. An amplifier and discriminator on each strip determines 
which strips were hit. For a uniform distribution of a variable x in an interval 
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(-a/2, a/2), 



1 f a/2 a 2 
(x z ) = - / (in 2 = — 

a -l-a/2 iZ 



so the nominal resolution orthogonal to the strip direction is about 80 fim/vl2. 
The resolution in practice is very similar. Measurement of the other coordinate 
is obtained by using small-angle stereo, i.e., by placing the strips of adjacent 
layers at a small angle. The innermost layers of the silicon detectors will 
use pixels of about 50 ^m x 300 ^m rather than strips. These are still under 
development. 

Central Tracker: Given the high multiplicity at the LHC, it is a difficult 
problem in pattern recognition to combine the right hits to find the real tracks. 
Having many layers in the tracking system helps, but it is prohibitive in terms 
of both cost and material to make many silicon layers. In ATLAS the outer 
portion of the tracker uses 70 layers of straw tubes. These have a thin con- 
ducting outer shell and a fine wire in the center, with a high voltage between 
them. Charged particles produce ionization in the gas, which is chosen to give 
a constant drift velocity for the produced electrons. These electrons drift to 
the central wire, where the high field produces an avalanche with a typical gas 
gain of <~ 10 4 . The time of this signal is measured, and the drift velocity is 
used to convert this to a drift distance with a resolution of order 100 ^m. 

The ATLAS tracker is a "Transition Radiation Tracker." It also incor- 
porates a transition radiation detector, which detects the X-rays emitted by 
a charged particle passing through a dielectric interface. This radiation is 
proportional to 7 and is useful for identifying low-py electrons. 

The CMS outer tracker uses gas microstrip detectors, which are still being 
developed. The ionization is produced in a layer of gas, but it is detected 
using closely spaced strips on a substrate read out in a manner similar to 
silicon strips. There are fewer layers than in the ATLAS tracker, but each 
layer has better position resolution and lower occupancy. 

The purpose of the central tracker is to determine the momenta of charged 
particles by measuring their curvature in the central solenoidal magnetic field. 
As a consequence of the Lorentz force, F = ev x B, a charged particle in a 
uniform solenoidal magnetic field B follows a helix with a radius of curvature 

R= PT 



0.3B ' 

with R in meters, B in Tesla, and pr in GeV. Simple geometry shows that the 
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resulting sagitta|| s for a radial length L -C R is given by 

L 2 _ 0.3BL 2 

S ~ 8R ~ 8p T 

with s and L in meters and pr in GeV. For pt — 100 GeV, B = 2T, and 
L = 1 m, the sagitta is s — 750 (im. The resolution depends on the chamber 
layout and position resolution and on the multiple scattering in the chamber 
material. ATLAS has B = 2 T, typical for most detectors, and R = 1 m, giving 

A ^~. 0.7 frfeU 0.014, 



Pt VlTeV, 

where the constant term comes from multiple scattering and is added in quad- 
rature. The resolution for tracks beyond the corner of the solenoid degrades 
like 1/ sin 2 6. CMS has a very high field, B = 4T, and also a larger radius, 
giving 

0.1 {f, 



Apr / PT 

p T ~ ' UTeV 



again degrading beyond the corner of the solenoid. ATLAS has a comparable 
resolution for muons only using its muon system. 

Electromagnetic Calorimeter: Precision electromagnetic calorimetry has 
been used in several e + e~ detectors, but ATLAS and CMS are the first hadron- 
collider detectors to have it. The demand for very high resolution is driven by 
the search for h — > 77; it will be seen in Section |^ that this gives a narrow 
peak on a large 77 continuum background. The ability to separate e and 7 
depends on the tracker, so the useful calorimeter coverage is also \t)\ < 2.5. 

Any calorimeter creates a shower in some dense material and uses the 
total charge or light output from this shower to determine the energy of the 
initiating particle. The energy is divided among more and more particles 
until it is completely absorbed. It is important to realize that electromagnetic 
interactions with a dense material like lead are much stronger than hadronic 
ones: electromagnetic interactions scale like Z 2 while hadronic ones scale like 
A 2 / 3 . For lead, the radiation length X , the distance in which a high energy 
electron loses all but 1/e of its energy, is 0.56 cm, while the inelastic hadronic 
interaction length A is 17.1cm. (For a light material like aluminum Xq = 
8.9 cm and A = 39.4 cm.) Thus, an electromagnetic calorimeter ~ 25Ao thick 
will contain almost all of the shower from high energy electrons or photons 
while absorbing little hadronic energy. Thus the electromagnetic calorimeter 
is always in front of the hadronic one. 

6 The sagitta, a standard term in elementary geometry, is the maximum separation be- 
tween the arc of a circle and its chord, the straight line between its ends. 
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The ATLAS electromagnetic calorimeter uses lead plates with gaps filled 
with liquid argon. Electrons from the shower drift under high voltage through 
the liquid argon and are collected on readout pads. The energy is proportional 
to the total charge, with an energy resolution 

AE 10% n M 
— w — = © 0.5% 

E y/E 

with E measured in GeV. The first term here comes from the shower multi- 
plicity and the fluctuations in sampling it: the multiplicity N of particles in 
the shower is proportional to the energy, and the Poisson fluctuation in N is 
1/VA- Note that this term gives the same resolution for one particle or from 
several with the same total energy, the errors being added in quadrature, as 
one would expect for a calorimeter. The small constant term arises from many 
sources and is added in quadrature. The ATLAS calorimeter can also use the 
position of the shower as a function of depth to measure the direction of a 
photon with an accuracy of about 50 mr/\/E. 

CMS uses a dense, transparent crystal, PbW04, both to create the shower 
and to convert it into scintillation light that can be detected by photodiodes. 
Because there are no inert lead plates, the whole shower can be measured, the 
sampling fluctuations are reduced, and the resolution is therefore better, 

¥->«■ 

with E measured in GeV. While crystals give better energy resolution, they 
make it harder to achieve fine segmentation and good pointing accuracy, and 
controlling the crystal quality is not trivial. At low luminosity, CMS will rely 
on tracking to determine the vertex and hence the photon direction. At high 
luminosity, it will add a preshower detector to measure the starting position 
of the shower and provide directional information at the cost of some energy 
resolution. 

Hadron Calorimeter: The hadron calorimeter follows the electromagnetic 
one and measures the energy of both charged and neutral hadrons. Hadronic 
showers have intrinsically larger fluctuations than electromagnetic ones, and 
additional errors are introduced for jets by the clustering algorithm. Both 
the ATLAS and the CMS central hadron calorimeters use steel plates (which 
also act as the magnetic flux return) and sample the hadronic showers with 
scintillators read out by photomultiplier tubes. (The ATLAS endcap hadron 
calorimeter uses copper plates and liquid argon.) The resulting jet resolutions 
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when combined with the electromagnetic calorimeters are roughly 
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Ve 



with E measured in GeV. The first term comes from the shower multiplicity 
and the fluctuations in sampling it, as for the electromagnetic calorimeter. The 
second term is added in quadrature and comes, e.g., from the fact that there 
are fluctuations in the fraction of the energy carried by 7r°'s and by charged 7r's, 
and the calorimeter responds differently to these.f] The forward calorimeters 
cover 3 < \rj\ < 5 with cruder energy resolution since E px in this region. 

Solenoid: The central trackers require solenoids to provide the magnetic 
field. The 2 Tesla solenoid in ATLAS is thin, < lAo, so it can be placed in 
front of the electromagnetic calorimeter without degrading its resolution too 
much. The 4 Tesla solenoid in CMS must be much thicker and so must be 
placed outside the hadron calorimeter. 

Muon System: Muons radiate much less than electrons and so penetrate 
the whole calorimeter with small energy losses, at least for energies below 
the TeV scale. ATLAS makes its precise muon measurement with an air-core 
toroidal magnet outside the calorimeter and three groups of tracking chambers 
to measure the resulting sagitta. The barrel toroid is made of eight lumped 
coils, which can be seen in Figure [| Lumped coils are needed to allow chambers 
to be placed within the toroid but give a rather complex, non-uniform field. 
The endcap toroid has one group of chambers in front of it and two groups 
behind it. The resolution pt for large pt in the central region is 



it remains quite good up to 77 = 2.5 because the endcap toroid gives a B oc 1/r 
at small radius, thus increasing the bending power at large 77. 

CMS relies on its central tracking to achieve comparable resolution for 
both muons and hadrons, making only a ~ 20% measurement in the exter- 
nal muon system, which utilizes the iron flux return of the central solenoid 
for its magnetic field. The CMS muon system provides triggering and muon 
identification rather than a precise measurement. 

Trigger: The trigger systems are crucial for the LHC experiments. The 
40 MHz interaction rate must be reduced by about a factor of 10 6 before the 

c It is possible to build calorimeters like that in the ZEUS detector at HERA with nearly 
equal electromagnetic and hadronic responses, but they generally have poorer electromag- 
netic resolution, the primary emphasis for both ATLAS and CMS. 
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Table 2: A set of possible Level 1 and Level 2 triggers for ATLASB at C = 10 34 cm _2 s _1 . 
The fix thresholds are not yet known but should be ~ 100 GeV; they will be set to give the 
indicated rates. 



Trigger Requirement 


LVL1 Rate 


LVL2 Rate 




(kHz) 


(kHz) 


> 1 muon, p T > 20 GeV 


4 




> 1 isolated fi, p T > 20 GeV 




0.2 


> 1 n,pr > 40 GeV 




0.1 


> 1 isolated e.m. cluster, Et > 30 GeV 


20 




> 1 electron E T > 30 GeV 




0.3 


> 1 isol. e.m. cluster, E T > 30 GeV 




0.1 


> 2 muons, px > 6 GeV 


1 




> 2 n, PT > 10 GeV 




0.1 


> 2 isolated fi, px > 6 GeV 




0.01 


> 2 isolated e.m. clusters, Et > 20 GeV 


4 




> 2 e or 7, p T > 20 GeV 




0.2 


> 1 jet, E T > 150 GeV 


3 




> 1 jet, E T > 300 GeV 




0.1 


> 3 jets, E T > 150 GeV 




0.04 


Missing energy ($r) 


1 


0.1 


Prescalcd triggers 


5 


0.1 


Total 


38 


1.4 



events can be saved to tape for future analysis, and this obviously limits the 
physics that can be studied. The trigger is divided into three levels. Level 1 is 
hardware-based, synchronous with the beam clock, and deadtimeless. The data 
from every detector element must be saved for about 2 /is until a Level 1 deci- 
sion can be made. This decision is based on fast sums of predefined clusters in 
the electromagnetic and hadronic calorimeter and on hits in roads correspond- 
ing to stiff tracks in the muon system. The thresholds on these are adjusted to 
reduce the rate to about 10 4 Hz. Level 2 refines the selection made at Level 1 
by using the full granularity of the detector and by combining measurements 
from more than one subsystem in "regions of interest" found by Level 1 trigger. 
This allows one, e.g., to determine the pt of an electron candidate more accu- 
rately both by using more detailed calorimeter information and by comparing 
it with tracking information. Finally, at Level 3 the whole detector is read out 
into a computer farm, which can run off-line analysis code and save events at 
roughly 100 Hz. A list of possible triggers for ATLAS and their Level- 1 and 
Level-2 rates is shown in Table @. 
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4-2 Measuring the Standard Model Quanta 

Any SUSY or other new particle will be produced at the LHC with px M, so 
its decay products will be widely distributed in phase space. It will either decay 
into quanta of the Standard Model — quark or gluon jets, charged leptons, 
neutrinos, or photons — or it will be stable and escape the detector like the 
lightest SUSY particle Xi if R parity is conserved. ATLAS and CMS are 
designed to detect such signals, in many cases with redundant measurements. 

Jets: Jets, the dominant signal at large pr, are measured as clusters of 
energy in the electromagnetic and hadronic calorimeters. There will also be 
multiple charged tracks connecting this cluster and the vertex. The jet energy 
resolution tends to be dominated more by uncertainties associated with jet 
clustering and QCD radiation than by detector performance. 

b Jets: The most important use of the vertex detector is to tag b jets. 
Most of the tracks in an event will point back to the primary vertex, but those 
from a B hadron will have a distance of closest approach characteristic of the 
B lifetime, i.e., ct as 465/im. Tagging is not easier for a highly relativistic 
B hadron: while the typical distance traveled by the B is jct, the typical 
opening angles of the decay tracks are 0(1/7), so the typical distance of closest 
approach to the primary vertex is independent of 7 provided j3 w 1. Calculating 
the tagging efficiency requires detailed simulation and obviously involves a 
tradeoff between efficiency and background rejection. For a_L% mistagging 



Photons: An isolated photon is identified as a cluster contained in the elec- 
tromagnetic calorimeter with a radius of order the radiation length Xq with 
no hadronic energy behind it and no high-p^ charged track near it. (Actually, 
the trackers in ATLAS and CMS both contain a significant fraction of IXq, so 
the probability that the photon converts to an e + e _ pair in the tracker is not 
small.) Isolation criteria are sufficient to reject jets by a factor of several thou- 
sand but still leave some background from those jets in which one or more 7r°'s 
carry most of the jet energy. An additional rejection can be obtained by using 
a preradiator to count the number of photons and/or by using detailed shower 
shape cuts. Thus, because of the very good ATLAS and CMS electromagnetic 
calorimeters, a jet rejection of ~ 10 4 can be achieved with a photon efficiency 
of order 90%. Hence, the background for h — > 77 should be dominated by the 
real QCD 77 continuum. Non-isolated photons of course cannot be separated 
from the much larger rate for 7r° — > 77 in jets. 

Electrons: An electron gives an electromagnetic shower like a photon but 
has a track with p as E pointing to it. Because of this extra constraint, isolated 
electrons with pr <; lOGeV and \r)\ < 2.5 can be identified with an efficiency 
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of order 90% with a jet rejection of ~ 10 5 , so that the background for isolated 
electrons is dominated in almost all cases by real electrons. Electrons within 
jets are much more difficult to identify. 

Taus: A r decaying into a lepton is difficult to distinguish from a prompt 
lcpton. A r decaying into hadrons can be identified as a narrow hadronic jet 
with one charged track or three tracks with a charge ±1 and a mass < M T . 
The background from QCD jets is significant, and detailed study is required 
to develop cuts appropriate for any particular case. 

Muons: A muon is identified by a charged track in the central tracker and 
a matching charged track in the external muon system. The energy deposition 
in the calorimeter is small (~ 2 GeV) for energies below about 1 TeV. At higher 
energies bremsstrahlung and e + e _ pair production become significant, and the 
energy deposited in the calorimeter needs to be considered. 

Missing Energy: In hadron colliders the total missing energy is completely 
dominated by the loss of low-px particles in the beam pipe, so only the missing 
transverse energy $t can be used to detect neutrinos or xT s - This is measured 
by summing all the calorimeters plus any observed muons. The resolution is 
dominated by non-Gaussian tails and cracks in the calorimeter; detailed studies 
indicate that real neutrinos dominate over the instrumental background at least 
for fJ T > 100 GeV. 

5 Inclusive SUSY Measurements at LHC 

If SUSY is indeed the right new physics at the electroweak scale, the first task 
of the LHC will be to detect a deviation from Standard Model predictions 
characteristic of SUSY. The ability to do so is clearly model dependent. For 
example, if all SUSY particles were nearly degenerate in mass, then they would 
decay into very soft jets or leptons plus an invisible x\, and nothing would be 
observable. Fortunately, such a degenerate spectrum does not occur in any 
reasonable model. 

5.1 Simulation of SUSY Signatures 

Most recent studies of SUSY signatures at the LHC have assumed the minimal 
supergravity (SUGRA) model.L3 The SUGRA model is a special case of the 
Minimal Supersymmetric Standard Model (MSSM), with two Higgs doublets 
and a SUSY partner for each Standard Model one, grand unification at some 
scale A/gut, and soft SUSY breaking terms added by hand assuming i?-parity 
conservation: 

-£ soft = A u h u QH u u c + A d h d QH d d c + A l h l LH d l c 
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-B^{H d H u + h.c.) + M 2 H \H d \ 2 



~> f\ i 71 ^2 I ~c 1 2 



+M^JH U \ Z + Mf\L\ J +M 
+Af||Q| 2 + Aif |u c | 2 + M||d c | 2 

+ Im x BB + X -M 2 WW + X -M^g. 

In SUGRA the soft breaking terms are assumed to be communicated from 
the SUSY breaking sector by gravity and so to be universal at Mqut- The 
resulting minimal set of parameters is: 

• 77io : the common SUSY- breaking mass of all squarks, sleptons, and Higgs 
bosons. 

• m i/2 : the common SUSY-breaking mass of all gauginos. 

• Aq: the common SUSY-breaking trilinear coupling. 

• £>/i: the SUSY-breaking bilinear coupling. 

• fi: the SUSY-conserving Higgsino mass. 

All of these parameters, including the SUSY-conserving parameter /j, should 
be of order the weak scale. A limitation of the SUGRA model is the absence 
of any understanding of why /i should be of order the weak scale or why R 
parity should be conserved. For a more detailed discussion, see the lectures by 
DawsonB in these Proceedings and references therein. 

The SUGRA model defines the SUSY breaking parameters at the GUT 
scale. All of these parameters are essentially couplings and so obey renormal- 
ization group equations (RGE's). The RGE's in the SUGRA model involve 
26 coupled partial differential equations, which have been studied by various 
authors;Ell an example is shown in Figure [l^. ISAJET implements a self- 
consistent solution of these RGE's between the weak and the GUT scale. The 
first step is to run a truncated set of six equations from Mz to the GUT scale 
where the gauge couplings g\ and g 2 meet using approximate SUSY mass scale: 

^1 - i f_3 _ 1\T,\„3 



dt 



16^ v 5 N f)9i + 2-loop terms 

^2 i (K _ ATAn3 

dt ^ 



dt 



(5 — Nf)g 2 + 2-loop terms 
Nf)g% + 2-loop terms 
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M = 300 GeV, M 1/2 =100 GeV, A = 
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Figure 13: Evolution of SUSY masses from the GUT to the electroweak scale. I 



Once Mqut is determined, the universal SUGRA boundary conditions are 
imposed, and the full set of 26 RGE's are run back to the weak scale using 
Runge-Kutta step-by-step integration so that mass thresholds can be properly 
taken into account. The Clebsch- Gordon coefficients in these equations are 
such that the Higgs mass is driven negative, breaking electroweak symmetry 
but not charge or color. The Higgs effective potential is determined, and the 
GUT scale parameters B and fi 2 are determined in terms of the weak scale 
parameters Mz and tan (3 — vi jv\ . The whole procedure is then iterated until 
a self-consistent solution is obtained. The final result is to express the masses 
of all 32 SUSY particles plus all the mixing parameters in terms of just four 
parameters plus sgn/i = ±1: 

• 77io : common scalar mass at Mqut- 

• common gaugino mass at Mqut- 

• Aq\ common trilinear coupling at A/qut- 

• tan/3 = vijv\\ Ratio of VEV's at M z . 

• sgn fi = ±1. 

In the SUGRA model the SUSY masses are mainly determined by mo and 
rrii/2, while tan/3 and sgn/i = ±1 mainly affect the Higgs sector. Aq is not 
very important for weak-scale physics; while A t , Ai,, and A T are important for 
third-generation sparticles, they turn out to be only weakly dependent on Aq 
over most of the parameter range. Hence it seems to be sufficient to scan the 
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too-TOi/2 plane for a few values of tan [3 and sgn/i = ±1 and one value of Aq, 
say A Q = 0. 

Since the same mo is used for all scalar particles, it must be that > 
so that charge and color are not broken. The SUGRA model is only possible 
because the top quark is heavy: it turns out that the large value of y t drives the 
Higgs mass-squared negative, breaking electroweak symmetry but not color or 
charge, as illustrated in Figure |l3|. SUGRA is surely not the final answer: it 
sheds no light on fermion masses, CP violation, etc. But it is a self-consistent 
framework representative of a large class of models, and it might even be close 
to the truth, so it seems worthy of serious study. Other models might be easier. 
In gauge-mediated models, the lightest SUSY particle is the gravitino; if the 
next lightest SUSY particle is the Xx> then Xx ~ * ^7 decays can be used to 
tag SUSY events with two hard photons. In i?-parity violating models, the 
Xi can decay either into three leptons or into three quarks; decays into both 
would lead to weak-scale proton decay. In the first case the leptons give a good 
signature. In the second, there presumably are still leptons fxom the cascade 
decays and it is possible to kinematically reconstruct masses.Ea 

The solution of the renormalization group equations for the SUGRA model 
is built into IS A JET. eJ The numerical solution of these equations uses Runge- 
Kutta step-by-step integration so that the thresholds corresponding to the 
various SUSY masses can be included in a self-consistent way. First, a trun- 
cated set of equations is used to determine a first estimate of the GUT scale, 
defined as the scale at which a% and ai meet. The GUT boundary conditions 
are then imposed, and the full set of equations is run back to the weak scale, 
freezing out mass parameters at their own scales. The 1-loop Higgs effective 
potential, including the SUSY masses, is computed, and the parameters B and 
/j, 2 are eliminated in favor of M§ and tan/3. Some optimization of the scale 
choice is made; this is equivalent to including some 2-loop contributions. The 
equations are then iterated until a self-consistent solution is found. Once the 
renormalization equations have been solved, the sfermion and gaugino mixing 
matrices are computed, and all the branching ratios for SUSY particles are 
computed. In earlier versions mixings in the b and f sectors were ignored, 
limiting the program to tan (3 < 10; this restriction has recently been removed. 
The branching ratio calculations use the correct matrix elements, but at present 
phase space is used in the actual event generation for technical convenience. 
Thus, for example, as a squark mass is varied from just below to just above 
the gluino mass, the branching ratio behaves sensibly but the event structure 
changes discontinuously. 

PYTHIA0 uses approximate formulas rather than solving the renormal- 
ization group equations, or it can take masses from an external calculation. It 
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Figure 14: Gluino and squark masseso in the minimal SUGRA model in the plane 
for tan/3 = 2, 10, sgn/i = ±1, and Ao = 0. The "bricked" regions are excluded theoretically, 
while the cross-hatched regions are excluded by experiment. The dots correspond to the five 
LHC points described in the text. 

treats the branching ratios in a similar way. 

Figures |lj and |l5| show contour plots of various SUSY masses in the mop; 
toi/2 plane for tan/3 = 2,10, sgn/x = ±1, and A = from ISAJET 7.22.1 
The cross-hatched regions are excluded by experiment. The bricked regions at 
small mo are excluded by the requirement that the Xi rather than the fi be the 
lightest SUSY particle. The bricked regions small miy 2 and tan/3 = 10 were 
excluded in ISAJET 7.22 by the absence of electroweak symmetry breaking. It 
turns out that the size of this excluded region is very sensitive to the scale at 
which the effective potential is minimized; changes in recent versions of ISAJET 
intended to make the results more stable for tan/3 ^S> 10 have significantly 
reduced this region. 



34 



800 



% 600 



J 400 
S 

200 


800 
« 600 

M 

J 400 
6 

200 





t 1 1 1" 1 1 1 I ■ 

; 4 


tin/3=2 - 

m<o : 






#,(700) — 
T E (1400) I 
#,(500) — 


j^ E (200) 


'; |T S (1000) 


' #,(300) ~ 












i #,(100) 












lllljlllt 

i ' b); 


tJm/3=2 1 

m>o : 






#,(700) — 
T E (1400) - 

#,(500) 


p • 


'; :U(1000) 
:T.(600) ; 


' #,(300) 


P 'T,<200) 

• 










i #,(100) - 
















500 1000 1500 2000 
m (GeV) 



500 1000 1500 2000 
m (GeV) 



Figure 15: Gaugino and slepton massestj in the minimal SUGRA model in the mo-ni]^ plane 
for tan/3 = 2, 10, sgn fi = ±1, and Aq = 0. The "bricked" regions are excluded theoretically, 
while the cross-hatched regions are excluded by experiment. The dots correspond to the five 
LHC points described in the text. 



Figures |l4j and |l5| illustrate a number of general features of the SUGRA 
mass spectrum: 



• Gluino and gaugino mass depend mainly on m^. 

• Slepton masses depend mainly on mg. 



Squark masses depend mainly on . ttVq + 4m? 



1/2- 



• Mq > 0.9Mg. 



M^o fa 0.5M^o ps 0.5M~± ~ 0.5 x 0.3M 5 . 

X-l X-2 Xl y 
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The last point is more general than the SUGRA model. It means that there 
is a large energy release at each step in the cascade decays of SUSY particles. 
If all the SUSY particles were nearly degenerate, they would be much more 
difficult to detect. 

5.2 Reach of SUSY Signatures 

Recall that a g or q is produced at the LHC with pt ~ M and decays into 
jets, possible leptons, and a Xii which is neutral and weakly interacting and so 
escapes the detector. Thus the most generic prediction of (i?-parity conserving 
weak scale) SUSY is an excess of events with multiple jets plus missing energy 
compared with the Standard Model sources, i.e., W, Z , and heavy quark pro- 
duction and mismeasurement of QCD jet events. The first task is to determine 
whether such an excess exists. 

The standard requirement for discovery of a new phenomenon is a signif- 
icance of at least 5a. That is, the probability that the background fluctuates 
up to the observed signal should be less than the tail of a Gaussian distribution 
beyond 5tr, i.e., 5.7 x 10 . This may seem overly conservative but is essential 
because one always looks at many different distributions with different cuts, 
and one of them is likely to have an unlikely fluctuation. For large numbers of 
events, the 5tr requirement is equivalent to 

S/VB > 5 

where S and B are the number of signal and background events respectively. 
For small numbers, Poisson probabilities should be used. Of course the S/B 
ratio, or more properly the error on determining what the background should 
be, must also be considered. Ruling out the existence of a signal is less de- 
manding, and limits-are generally quoted for 90% or 95% confidence. 

The approachEm for determining the LHC reach in the SUGRA parameter 
space is to scan the mo-m^ plane for selected values of the other parameters, 
generating a sample of SUSY events for each choice of parameters. These 
consist mainly of g and q production, but all processes are included. Since 
one event typically takes about 0.5 s on an HP-735, the feasible data samples 
correspond to ~ lOfb -1 for detailed studies but much less for such a scan. 
It is clearly not possible to generate a representative sample of the Standard 
Model total cross section. Instead, high-p^ events which potentially can give 
large fx, namely 

• W+jets, W -> tv,rv 

• Z + jets, Z — > vD, tt. 
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• tt. 

• QCD jets, including g — > 66, cc branching and decay. 

These samples are generated using several approximately equal intervals of 
logpT for the primary hard scattering. The event generator of course pro- 
duces not just the hard scattering but also parton showers, hadronization of 
the partons into jets, and beam jets. The studies described here generally 
have assumed low luminosity, C = 10 33 cm - 2 s -1 , and have therefore neglected 
pilcup from overlapping events. 

All events are passed through a toy detector simulation. This takes into ac- 
count the overall coverage and Gaussian resolutions but not cracks, resolution 
tails, multiple scattering, or many other effects. It is possible to take all these 
effects into account, but the detector simulation then requires more than 1 hr 
per event. The toy simulation is not adequate to determine the background 
from mismeasured QCD jets. More detailed studies show that this background 
is less than that from real Standard Model neutrinos for Ifix ft 100 GeV, and 
this cut will generally be made whenever fx is used. 

Jets are found using a simple fixed-cone algorithm. That is, the highest 
remaining unused cell of the calorimeter is found, and the total Et is summed 
in a cone in 

r = ^(A V r + (W 2 , 

which is equivalent to a polar angle at r\ = but is z-boost invariant and 
so behaves properly in the forward direction. Generally R = 0.4 is used for 
complex multijet events as a compromise between identifying nearby jets and 
containing all of the jet energy. Electrons and muons are treated equivalently, 
requiring isolation Et < 5 GeV in cone R = 0.3 both to suppress the b, c — » £X 
background and to permit e identification by track/shower matching. The lat- 
ter requires a more detailed simulation to implement, so generator information 
is used for lepton identification. 

Figure |l6| summarizes the reach of the LHC to observe SUSY in various 
channels in the mo — mi/2 plane for two representative values of tan /3, sgn /j, = 
±1, and Ao = 0. The reach limits are based on a 5a signal after 10/6 -1 , 
corresponding to one year at low luminosity. 

0£, l£: These curves show the reach in the basic channel, multiple jets plus 
missing energy. The 0£ curve includes a veto on muons and isolated electrons, 
while the 1£ curve requires a lepton. The lepton veto improves the S/B ratio 
for low masses, but at high masses so many leptons are produced that requiring 
a lepton improves the reach. For both sets of curves the following cuts are made 
to reject the Standard Model background and to enhance the acceptance for 
heavy particles produced with pt ~ M : 
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m (GeV) m (GeV) 

Figure 16: SUGRA reach in the mj-m]^ plane for tan j3 = 2,10, sgn/x = ±1, and Aq = 
for 10 fb" 1 luminosity at the LHC.B The "bricked" regions are excluded theoretically, while 
the cross-hatched regions are excluded by experiment. The dots correspond to the points 
selected by the LHCC. 



• Njcts > 2 with E T > 100 GeV, \r)\ < 3. 

• Closest jet j c has 30° < A(j>(fi T ,j c ) < 90°. 

• S T > 0.2 

• f T >E c T ,E T { ]l ),E T {j 2 )>E c T . 

Here Ej. is a cut which is adjusted at each point to optimize S/y/~B, and the 
optimum value is used to define the reach. The variable St is the transverse 

38 



sphericity or circularity, which is defined as 



St = 



2X 2 



Ai + A 2 



where Ai > A 2 are the eigenvalues of the transverse sphericity tensor 



This cut selects "round" events characteristic of heavy particle production, but 
it is highly correlated with the multijet and other cuts. 

SS: This curve shows the reach in the like-sign dilcpton channel, £ ± £ ± , 
I = e,fi. The leptons are required to have pt > 20 GeV and \r]\ < 2.5, and to 
satisfy the isolation criterion Et < 5 GeV in a cone R = 0.3. Since the gluino 
is a Majorana fermion, it is its own antiparticle and so has equal branching 
ratios into £ + X and £~X, e.g., through xt cascade decays. (There may be 
other SUSY sources of like-sign dileptons.) The dominant Standard Model 
isolated dilepton backgrounds, tt, Drell-Yan, and W + W~ , only give opposite- 
sign dileptons. There are Standard Model like-sign dilepton backgrounds, e.g., 
from tt production with t — > C + X and t — > b — > £ + X, but these will normally 
fail the isolation test. 

OS: SUSY can also give opposite-sign dileptons £ + £~, £ = e,/i. e.g., 
from X2 ~~ > X?^ + ^~- The same lepton cuts are made as for the SS curves, 
and a Z mass cut is also made for identical flavor leptons. Opposite-sign 
dilcpton decays are enhanced at low mo, for which the sleptons are light and the 
decay can proceed via X2 — > ^ — > Xi^ + ^~ with substantial branching ratios. 
While the Standard Model backgrounds are larger in this channel than for SS, 
the statistical reach is comparable; the OS channel also provides important 
independent information. 

3£: SUSY can produce trilepton events from a variety of sources, including 
the decay of one gluino or squark via xt — > x\^ ±v an d the other through 
X2 ~~ * Xi^ + <?~. The Standard Model background for three isolated leptons is 
fairly small, so the reach in this channel is comparable to the dilepton channels 
even though the total branching ratio is smaller. 

31, Oj, 2£,0j: These channels require two or three leptons with the same 
cuts as before, $t, arid no jets with pt > 25 GeV and \rj\ < 3. The jet veto 
is designed to select the direct production of gaugino or slepton pairs. These 
channels are the best way of searching for SUSY at the Tevatron, where the lim- 
ited energy suppresses the production of the heavier gluinos and squarks. The 
search range is limited by competition from X2 ~ * Xi^, Xi^ 1 once mi/2 is large 




n 
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enough that these are kinematically allowed.u Even at smaller mxii the branch- 
ing ratio can be suppressed by interference between the virtual Z and slepton 
exchange graphs, leading to the holes in the reach seen in Figure |1| Neverthe- 
less, these channels would provide useful additional information should they 
be observed. 

By comparing Figure |l^ and the mass contours in Figure [li], one can 
see that the LHC can search the whole SUGRA parameter space at the So- 
level for gluino and squark masses up to about 2TeV with only 10fb _1 rif 
luminosity. Similar conclusions have been found by the ATLASEil and CMSEj 
Collaborations using the more general MSSM model and a more realistic pa- 
rameterization of the detectors. In addition, various multilepton signatures 
can be observed for gluino and squark masses up to about 1 TeV, i.e., over 
the whole range favored by fine-tuning arguments. The multilepton signatures 
with a jet veto are limited to relatively small values of mi/2 or mo- Thus, if 
SUSY exists at the weak scale, ATLAS and CMS should observe characteristic 
deviations from the Standard Model after one year of operation at only 10% 
of design luminosity. The LHC should either find SUSY or exclude it. Only 
experiment will decide whether SUSY at the weak scale is a crucial element in 
physics or an interesting exercise in mathematics. 

5.3 Introduction to Precision Measurements 

While observing signatures characteristic of SUSY at the LHC would be one of 
the most exciting developments in particle physics of all time, it is important 
to be able to determine the masses and other parameters of SUSY particles 
and thus to get a handle on the underlying dynamics. If R parity is conserved, 
however, then every SUSY event is missing two Xi's, an d there are not enough 
kinematic constraints to determine their momenta. 

It may be useful to compare the SUSY case with the production of tt at 
the Tevatron: 

q + q -> t + t -> e + vb + q'q"b . 

In these events there is one missing v and hence three unknown kinematic vari- 
ables p u . To determine these, there are two measured components of $t and 
two additional constraints expressed as quadratic equations in the components 
of p y : 

(Pe+Pu) 2 = 
(Pe +Pu+ Pb) 2 = (Pq> + Pq" + Plf ■ 

Thus there is one more constraint than unknown; this is known in ancient 
bubble chamber terminology as a 1C fit. Using all the constraints, one can 
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Table 3: Parameters for the five LHC SUGRA points. 



Point 


m 
(GeV) 


%/2 

(GeV) 


A 
(GeV) 


tan/3 


sgn/x 


1 


400 


400 





2.0 


+ 


2 


400 


400 





10.0 


+ 


3 


200 


100 





2.0 




4 


800 


200 





10.0 


+ 


5 


100 


300 


300 


2.1 


+ 



fully reconstruct the event despite the missing neutrino. If top were produced 
singly, then one would have only one quadratic constraint, a 0C fit; in this 
case p v could still be reconstructed, but there would be a 2-fold ambiguity. Of 
course for ti production one can also reconstruct the 3-jet mass directly. 

For SUSY, there are two missing Xi s an d so six unknown momentum com- 
ponents in addition to the Xi mass. The SUSY signal contains many different 
processes; there is no simple constraints like the W mass in the top case, so 
there are only two constraints on the six variables from the two components of 
fx- Hence it is not possible to reconstruct the events in general. It is possible, 
however, to use the kinematic endpoints of various distributions to make a pre- 
cise determination of combinations of masses. This is simplest in a e + e~ (or 
/j + /i~) collider, where the SUSY particles are produced with a known beam 
energy, giving extra constraints, but it is also possible at the LHC, as will be 
seen in Section |^ below. 

While the possibility of making such precision measurements is quite gen- 
eral, which ones can be made depends on the assumed masses and branching 
ratios and so can be determined only by simulating in detail events for specific 
choices of the SUSY parameters. The LHC Committee (LHCC), the CERN 
committee overseeing the LHC experiments, selected the five SUGRA points 
listed in Table || for detailed study by the ATLAS and CMS collaborations. 
Point 3 is the "comparison point," selected so that every existing or proposed 
accelerator could discover something. At this point, the Tevatron would dis- 
cover winos and zinos, the LHC would discover gluinos and squarks, the NLC 
would discover sleptons, and LEP would have recently announced the discov- 
ery of a light Higgs boson with a mass of 68 GeV. Points 1 and 2 have gluino 
and squark masses of about 1 TeV and so test the reach of the LHC for such 
masses. Point 4 has large tooi so that sleptons and squarks are much heavier 
than gauginos and gluinos. It was also close to the boundary of the allowed 
electroweak symmetry breaking region with ISAJET 7.22, so that y, was quite 
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Table 4: Masses of the superpartners, in GeV, at the five LHC SUGRA points from 
ISAJET 7.22. The first and second generation squarks and sleptons are degenerate. 



Point 


1 


2 


3 


4 


5 


g 


1004 


1009 


298 


582 


767 


xt 


325 


321 


96 


147 


232 


~± 

X2 


764 


537 


272 


315 


518 


xl 

A. 1 


168 


168 


45 


80 


122 


xl 


326 


321 


97 


148 


233 


x ( l 


750 


519 


257 


290 


497 


xi 

/V-i 


766 


538 


273 


315 


521 


ul 


957 


963 


317 


918 


687 


UR 


925 


933 


313 


910 


664 


d L 


959 


966 


323 


921 


690 


d R 


921 


930 


314 


910 


662 


h 


643 


710 


264 


594 


489 


h 


924 


933 


329 


805 


717 


h 


854 


871 


278 


774 


633 


h 


922 


930 


314 


903 


663 


e-L 


490 


491 


216 


814 


239 




430 


431 


207 


805 


157 




486 


485 


207 


810 


230 


n 


430 


425 


206 


797 


157 


T-2 


490 


491 


216 


811 


239 




486 


483 


207 


806 


230 


h° 


111 


125 


68 


117 


104 


H° 


1046 


737 


379 


858 


638 


A 


1044 


737 


371 


859 


634 


H± 


1046 


741 


378 


862 


638 



small and there was large mixing between the gauginos and higgsinos. More 
recent versions of ISAJET find that Point 4 further from this boundary, so that 
ji is larger and the mixing of gauginos and higgsinos is smaller. Finally, Point 5 
was chosen to be in the center of the region giving the right amount of cold 
dark matter for cosmology.cHl Heavy stable xj's tend to overdose the universe; 
getting the right amount of cold dark matter generally requires enhancing the 
Xi annihilation cross section and hence having relatively light sleptons. 

The masses of the SUSY particles for these five points as calculated with 
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LHC Point 1 




M eff (GeV) 



Figure 17: LHC Point 1 signal and Standard Model backgrounds. El Open circles: SUSY 
signal. Solid circles: it. Triangles: W — » iv, TV. Downward triangles: Z — > i/p, tt. Squares: 
QCD jets. Histogram: sum of all backgrounds. 



ISAJET 7.22 are listed in Table ^. These masses and the corresponding branch- 
ing ratios are used in all the analyses described below. 

5.4 Effective Mass Analysis 

The SUSY reach limits discussed in Section || are based on just counting the 
number of events with some specified set of cuts. Because QCD corrections to 
hadronic cross sections are large, the signal expected in the Standard Model 
is somewhat uncertain. It is therefore desirable to measure for some variable 
a distribution which agrees with the Standard Model in some range and then 
deviates from it, thus giving a more convincing signal and also providing an 
estimate of the SUSY mass scale. 

What properties should such a variable have? At least in the SUGRA 
model, the squarks are never much lighter than the gluino, so gluino produc- 
tion, which is enhanced by color and spin factors, is always important. If the 
squarks are heavier than the gluino, then the dominant gluino decays will be 
g — > xqq, giving a minimum of four jets plus missing energy. If the gluino is 
heavier, then the decay chain g — > qq, q — > \q will dominate. In either case 
there will be at least four jets plus missing energy, more if one or more of the 
gauginos in the process decay hadronically. QCD cross sections fall rapidly 
with momentum transfer — the jet cross section at y/s = 14TeV falls over the 
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M eff (GeV) 



Figure 18: LHC Point 2 signal and Standard Model backgrounds. See Figure [IT] for defi- 
nitions of the symbols. 




M ell (GeV) 



Figure 19: LHC Point 3 signal and Standard Model backgrounds .Cj See Figure [17] for defi- 
nitions of the symbols. 



relevant range of px roughly like 

da 
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LHC Point 4 




M eff (GeV) 



Figure 20: LHC Point 4 signal and Standard Model backgrounds .cil See Figure [17] for defi- 
nitions of the symbols. 




M elf (GeV) 



Figure 21: LHC Point 5 signal and Standard Model backgrounds .r H See Figure [IT] for defi- 
nitions of the symbols. 

— so it is clearly important to compare signal and background at comparable 
Pt scales. The invariant mass of the produced system is not the best measure 
of this because it is too much influenced by possible soft jets at large rapidity. 
The scalar sum of the px's of the four hardest jets and the $t works well and 



45 



Table 5: The value of M e ff for which S = B compared to MsuSYi the lighter of the gluino 
and squark (ur) masses. Note that Point 3 is strongly influenced by the and jet px cuts. 



LHC Point Mgg (GeV) MgugY (GeV) Ratio 

1 1360 926 L4T" 

2 1420 928 1.53 



3 470 300 1.58 

4 980 586 1.67 
5 980 663 1.48 

will be called the "effective mass"0 

M eS = p TA + p Tj2 + PT.3 + Pt,4 + $t ■ 

Backgrounds from QCD processes with multiple jets and neutrinos from heavy 
quarks generally have missing energy small compared to the Q 2 scale of the 
event. To avoid such backgrounds the fx cut is made proportional to M c g, 

$ T > 0.2M cS , 

where the coefficient 0.2 was chosen after studying the SUSY and Standard 
Model Monte Carlo distributions. 

Several additional cuts were made for technical reasons: a missing energy 
cut IpT > 100 GeV to ensure that the Standard Model $t background is 
dominated by neutrinos rather than mismeasured jets; a jet cut p-r,j > 50 GeV 
to ensure that the jets were well identified and measured, and a cut pt,i > 
100 GeV on the hardest jet to limit the range of QCD background that had to 
be generated. These cuts require M c g > 350 GeV and so limit the sensitivity to 
very light SUSY particles. In addition, there is a cut on transverse sphericity 
St > 0.2 to select "round" events characteristic of SUSY, although this is 
highly correlated with the previous cuts. Finally, there is a veto on muons or 
isolated electrons with px > 15 GeV; this minimizes the background for SUSY 
masses comparable to the top mass but reduces the sensitivity for high masses. 

After all these cuts, the Standard Model background dominates the M e g 
distribution for low M e g, but the SUSY signal dominates by a factor of 5-10 
for large M e g for all of the LHC points except the comparison point, Point 3, 



as can be seen from Figures |17|-|21|. For Point 3, Figure 19, the SUSY signal 
is larger than the Standard Model background for all values of M e g allowed 
by the technical cuts described above. Observing such a change in shape from 
a curve dominated by Standard Model physics to one a factor of 5-10 larger 
would be convincing evidence for new physics. 
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Figure 22: Scatterplot of MsuSY = min(Mg,Ma) vs. M e ff for randomh^-phosen SUGRA 
models having the same light Higgs mass within ±3 GeV as LHC Point 5-r 1 ! 



The signal curves in Figures |17|-|21| clearly shift with the SUSY masses. 
Since the SUSY cross section is dominated by g and q production, it is natural 
to use 

Msusy = min(M 5 ,M SR ) 

as a measure of the mass scale. Table pi shows the points at which the signal 
and background points cross Figures [l7|— |2l|. Clearly these points scale quite 
well with Msusy ■ 

To see whether the scaling in Table [s] might be accidental, a comparison 
of M e ff was madeO between 100 random SUGRA models and Point 5. The 
models were generated with parameters uniformly distributed in the intervals 
100 < m < 500 GeV, 100 < m 1/2 < 500 GeV, -500 < A a < 500 GeV, 
1.8 < tan/3 < 12, and sgn/i = ±1. All 100 models were selected to have 
M/j within an assumed theoretical uncertainty of ±3 GeV from the Mh mass 
for Point 5. The value of Af c ff for each model was determined not by the 
intersection with the Standard Model background but by the peak of the signal, 
which is somewhat higher. (It is not at all obvious that this is the optimal 
procedure, but it is what has been done.) The scatter plot of the peak M e ff 
vs. Msusy for each model is shown in Figure E2, and the projection is shown 



in Figure 23. Evidently the scaling of M c g vs. Msusy works remarkably well 
for this random selection as well as for the five LHC points. While the scaling 
is physically plausible, it is not known how well it works for arbitrary SUSY 
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models. 

6 Precision Measurements with Exclusive Final States 

While M ef j seems to work quite well as a measure of the SUSY mass scale, 
it clearly averages over many final states and branching ratios, so it can only 
be a rough approximation. To do better, one needs to reconstruct specific 
final states. If R parity is conserved, then every SUSY event is missing two 
Xi's, so no masses can be reconstructed directly. It is possible, however, to 
determine precisely combinations of masses by finding endpoints of kinematic 
distributions in specific final states, starting at the bottom of the decay chains 
for the SUSY particles and working up.Elc3 For simple SUSY models such a 
SUGRA with only a few parameters, this approach can determine the model 
parameters with good accuracy, at least in favorable cases. Even for more 
complicated models it is a good starting point. This Section describes a number 
of such precision measurementscirEa for the five LHC SUGRA points. 

6.1 Measurement of Mix") - M(xT) 

First consider Point 3. This point has unusual branching ratios because M^ < 
Mg but Mq > Mr, so that g — > bb is very much enhanced: 

B(g -> y> + h.c.) = 89%, 
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LHC Point 3 




Figure 2ir—P, Jr l mass for Point 3 (solid) and Standard Model background (shaded, nearly 
invisible) CJ 



B(h -» x° 2 b) = 86% , 
B(x°2^x1t + n =33%. 

While these branching ratios are not typical, it is common for heavy flavors in 
SUSY decays to be comparable to light ones or even enhanced. 

The SUSY particles at this point are relatively light and so give fx in the 
range for which detector effects are not negligible. Hence, fx is not used in 
the event selection at this point. Instead, events are selected by requiring 

• e+£- pair with p T j > lOGeV, 77 < 2.5. 

• > 2 jets tagged as b quarks with pr > 15 GeV and r\ < 2. 

• No fJ T cut. 

making use of the large X2 ~~ > Xi^ + ^~ branching ratio. The details of the 
selection are certainly specific to this point, but it should be possible in general 
to use leptonic modes to observe SUSY particles in this low mass range. 

SUSY signal and Standard Model background events were simulated as 
described above, and events were selected with the criteria just listed. Then 
the £ + £~ mass distribution was plotted, including a 60% tagging efficiency 
for 6's and a 90% efficiency for electrons and muons. This mass distribution 
shows a spectacular edge at the M(x°) — M(x1) endpoint, Figure [24|. This 
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Figure 25: Mg+g- distribution at Point 4 for opposite-sign, same-flavor dileptons 
(solid), orrmosite-sign, opposite-flavor dileptons (dashed) and Standard Model background 



distribution reflects the strong signal production, the large branching ratios, 
and the distinctive signature, resulting in almost no Standard Model back- 
ground. The signal has huge statistics, and measuring the position of the edge 
is clearly easier than measuring the W mass at the Tevatron, since only the 
lepton resolution and not the global $t resolution enters. Since the latter has 
already achieved an error of about 40 MeV, it seems conservative to estimate 
that one could measure the position of this edge and determine MijQ) — M(x?) 
to - 50 MeV. 

The huge statistics from the low masses and the extra handle of the b tag 
make this measurement unusually easy. But there is a similar edge at Point 4 
plus a Z peak coming from decays of the heavier gauginos. The dominant 
background is from two independent Xi decays in SUSY events. This and the 
Standard Model backgrounds, e.g., from top decays, contribute equally to all 
combinations of leptons and so vanishes up to statistical fluctuations in the 
combination e + e~ + — e + /i~ — e~/i + . Figure ^5] shows a plot of this 

difference for Point 4. The fluctuations in the background reflect the limited 
Monte Carlo statistics. Given the number of signal events, one could measure 
A { M {xD ~ M i.Xi)) = ±1 GeV with lOftT 1 at this point. 

This method should work for any point for which the direct 
branching ratio is not too small. This is generally the case for M(x°) ~ 
200 GeV so that x° ~h Xi^: Xih- -But there is a region of m for sgn^i = 
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+1 where the interference of the Z and slepton exchange graphs makes the 
branching ratio small. This is the same effect that produces the holes in the 
3£, Oj reach curves in Figure 16. 



6.2 Reconstruction of g and bi 

The X2 ~* Xi^ + ^~ decays at Point 3 can be combined with 6 jets to reconstruct 
the g — > bib, b\ — > x\° decay chain. The trick is to select an £ + £~ pair near 
endpoint of the distribution. Then the xi must be soft in the £ + £~ rest frame, 
so that 



where M^o must be determined from the model. (A first approximation is 
M(x°) = 2M(x?), so that the £+£~ endpoint directly gives M(x?)-) The 
following cuts are made for this analysis: 

• > 2 jets tagged as b jets with p T > 15 GeV, rj < 2; 

• £+£- pair with 45 < M{££) < 55 GeV. 

Again, no use is made of $t- Then the inferred x.2 momentum is combined 
with one b jet to make M(pi), and a second b jet is added to make M(g). A 
scatter plot including all the possible combinations for each event is shown in 



Figure 26. The projections of this scatter plot are shown in Figure |27|. 

Figure |26| was made assuming the correct Xi mass, which is not directly 
measurable but must be determined from a global fit as described in Section |?] 
below. If M(xi) is varied about its nominal value, then the value of M(bi) 
extracted from Figure ^7] shifts linearly, while the value of M(g) — M(bi) 
remains constant. Thus the errors on these masses are estimated to be 

AM(6i) = ±1.5AM(x?) ± 3 GeV , 
A(M(g) - M(h)) = ±2 GeV . 

The fact that the difference A(M(<?) — M{b\)) is insensitive to the assumed 
mass is a simple consequence the kinematics and the fact that one b jet is 
soft. Exactly the same effect is familiar from D* — > Dir decays. In addition 
to the uncertainty from the x.i mass, we must worry about the calibration of 
the hadronic energy scale. This can be studied using Standard Model events. 
For example, one could use Z — > e + e~ events to calibrate the electromagnetic 
calorimeter and then use px balance in 7 + jets and Z + jets events to transfer 
this calibration to the hadronic calorimeter. This is an important issue for the 
detector collaborations but need not concern us here. 
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Figure 26: Scatter plot of M(b + \2) vs. M(bbx2) — M(bx2) assuming M-o known 
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Figure 27: Projections of Figure The dashed curve shows the mass difference for events 

in the b mass peak. 



6.3 Reconstruction of h — > bb 

For Point 5, the decay x° Xih nas a large branching ratio. A set of cuts 
similar to those made in the M e s analysis is used to select the SUSY events 
and reject the Standard Model background: 
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Figure 28: distribution (e* Point 5 (solid), a Gaussian fit to the peak, and the Standard 
Model background (shaded). r 1 ! 



• > 4 jets with p T > 50 GeV, p T ,i > 100 GeV; 

• Transverse sphericity St > 0.2; 

• Afcff > 800 GeV; 

• $ T > max(100GeV,0.2M off ). 

After these basic cuts, at least two jets are required to be tagged as Vs with 
the vertex detector and to have px,b > 25 GeV and r/i, < 2. The distribution 
of masses of all pairs of b jets is shown in Figure |2^. The width of the peak 
at the Higgs mass is set by clustering effects and resolution. The dominant 
background is other b jets from SUSY events, not Standard Model background. 
This peak is much easier to detect than h — ► 77 and would be the discovery 
mode of the Higgs boson at this point. 



With the clustering algorithm described in Section 5.2 the Higgs peak 
actually came out several GeV below the Higgs mass. This is due partly 
to missing neutrinos and partly to leakage of energy out of the jet cone. A 
correction, 

2.9 GeV) , 

In the 

real experiment, one would need to use the precision mass from h — > 77 to 
determine this correction. 



1.08 (E° hscivcd 



was derived from the Monte Carlo data and applied to make Figure 28 



53 



LHC Point 1 LHCPoint2 




The branching ratio for Xih is generally substantial if it is kinemat- 

ically allowed, i.e., for m 1 / 2 ^ 250 GeV. There are similar Higgs mass peaks 
for the two other LHC points in this range; see Figure ^9|. One can hope to 
observe either \2 ~^ Xv 1 or X2 ~* Xi^ + ^~ over most of the parameter space. 

It is also possible to reconstruct W — ► qq at Point 5. Since there is no b 
tag to select the relevant jets, the Standard Model backgrounds are large, and 
very hard cuts are needed.Ej 

6.4 Reconstruction of g + g — > q^q + q^q 

Just as x° — * Xi^ + ^~ was use d as the starting point to reconstruct partially 
61 and g at Point 3, one can use h — ► bb as the starting point for more com- 
plex analyses. For example, at Point 5, the gluino is somewhat heavier than 
the squarks but has a larger cross section because of color and spin factors, 
producing as one possible signal 

9 + 9^ q~Lq + qaq , 
9l -> xll Xihq, qR -> Xil ■ 

The quark jets from the q decays are hard, while the other jets are softer. These 
events can be selected by requiring in addition to the two b jets two and only 
two additional jets with px > 75 GeV. Then one of the two qbb combinations 
should come from a q decay, so the smaller of the two qbb masses should have 
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Figure 3&r-i Smaller of two bbj masses at Point 5 (solid) and Standard Model background 
(shaded) O The background statistics are limited by the Monte Carlo event sample. 

an endpoint at the kinematic limit for q — > y^q —> x^hq. This limit is readily 
found from two-body kinematics to be 

(M,7 X ) 2 = Ml + (M| - M|) 

M| + Ml - M| + J(Aq - Ml - Mp - 4Mp| ' 

X ^5 . 

The average of the ul and <Il masses gives M™ x = 506 GeV, which is consis- 
tent with the edge seen in Figure The estimated error on this combination 
of masses is about one bin or 40 GeV for 10 fb _1 . It is statistics limited and so 
should improve with more luminosity. 

6.5 £ + £~ Distribution at Point 5 

The £ + £~ mass distribution has already been used to measure M(xV) ~ M(xi) 
at Points 3 and 4. Consider now this mass distribution for Point 5, one of the 
points for which a x§ Xih signal has been observed. The event selection 
combines the same basic cuts as before plus the requirement of a lepton pair: 

• M cS > 800 GeV; 

• $ T > 0.2M eff ; 
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Figure 3trn M(£ + £ ) distribution for Point 5 (solid) and Standard Model background 
(shaded) IB 

• > 1 jet with p T > 100 GeV; 

• Isolated £ + £~ pair with p T > 10 GeV, r\ < 2; 

• Transverse sphericity St > 0.2. 

The opposite-sign, same-flavor dilepton distribution after these cuts is shown 
in Figure ||l] and obviously has a dramatic edge with very little background. 

It might at first seem that one should interpret this edge as a measure of 
M(x<j) — M(xi) from the kinematic limit of x\ ~ * Xi^ + £~ ■ But this interpre- 
tation is very implausible given the observation of h — ► bb peak at this point, 
since this would require a 3-body decay to compete with a 2-body one, whereas 
3-body phase space is much smaller than 2-body phase space. A much more 
plausible explanation is that another two-body decay, Xz ~ ¥ ^> gi yes r i se to 
the dilepton endpoint. In SUGRA one generally has M(£r) < M(£l), and in 
fact the observed endpoint in this case comes from 



Simple two-body kinematics for the X2 ~ ¥ an d @R ~ * Xi^ decays gives for 
the M{£ + £~) endpoint 
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Figure 32: Same as Figure Bll but fop-Point 5 (solid) and tuq = 120 GeV (dashed), and 
Standard Model background (shaded) Jfl| 



The estimated error on this combination of masses is about 1 GeV for 10 fb~ . 
It is statistics limited and so should improve with more luminosity. 

This example illustrates that the interpretation of edges of kinematic dis- 
tributions at the LHC may be ambiguous. If a signal for %° ~* Xih, h — > bb, 
had not already been observed at this point, it would not be clear whether 
the edge measured M(\2) — M(x1) from direct decays or the above combina- 
tion of three masses from decays through sleptons. Thus it is important to 
measure as many combinations of masses and distributions as possible so as 
to overconstrain models. 

To get the maximum information from the dilepton distribution at this 
point, one should fit at least the M(££), pr(££), pr(h), and pt(£2)/pt(£i) dis- 
tributions varying M(x2), M(£r), M(xi), and the pr distribution of This 
is certainly what one would do if the distribution in Figure [n] were real exper- 
imental data, but it requires more effort than seems warranted now. Instead, 
one additional sample of events was generated with mo = 120 GeV rather than 
100 GeV; this change increases M(£r) by 13 GeV while having small effects 
on all the other relevant masses. As a result there is a 2 GeV change in the 
location of the edge, as seen in Figure |32|. The changes in the pt distributions 
are small. The most significant change is in the distribution of the variable 
Pt(£2)/pt(£i), which is shown in Figure This makes physical sense: in- 
creasing the slepton mass reduces the phase space and hence the pr for the 
first lepton, while increasing the phase space and the pr for the second lepton 
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Figure 33: Ratio Pt{^2)/pt^i) for Point 5 (solid), mo = 120 GeV (dashed), and Standard 
Model background (shaded) .Cj 

from the slepton. It is clear that a change of this magnitude is easily detected; 
the estimated error on the slepton mass from such measurement is ~ 5 GeV on 
mo and ~ 3 GeV on the slepton mass. Changing the slepton mass also changes 
the \2 ~ * branching ratio, providing additional information. 

6.6 Measurement of M(g) - M(x°) and M(g) - M(xi) at Point 4 

Gluino production dominates at Point 4 since mo mi/2, so that the squarks 
are much heavier than the gluino. An analysis described previously for this 
point found a x§ — > x\t + i~ edge and hence a measure of M(x°) — ^(x?)- The 
jet multiplicity from gluino decay is large, so there is a lot of combinatorial 
background, pand reconstructing a gluino signal is not trivial. The strategy for 
this analysisca is to select 

9 + 9 xlm + xfqq 

using leptonic decays of both gauginos to identify them and so to reduce the 
combinatorial background. Then the right combination of jet-jet masses has a 
common endpoint since M(%°) ~ M(xt)- 

The event selection for this analysis imposes the following set of cuts: 

• 3 isolated leptons with pt > 20, 10, 10 GeV and \r)\ < 2.5, with one OS, 
SF lepton pair. 
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Figure 34: Jet-jet mass distribution for Point 4 with three leptons and four jets (solid) and 
for correct pairing (dashed) .E3 




• M < 72 GeV, the edge in Figure 25, for the OS, SF pair, consistent with 
X20 decay and not with possible Z backgrounds. 

• 4 jets with p T > 150, 120, 70, 40 GeV and \rj\ < 3.2. 
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• No additional jets with px > 40 GeV and \rj\ < 5 to minimize combina- 
torial backgrounds. 

• No $ T cut. 

After these cuts there 250 signal events with 30 gq background and 18 other 
SUSY and Standard Model background for 30 fb -1 of integrated luminosity. 
Thus the signal purity is quite good, and the signal efficiency, ~ 35%, is also 
good. 

There are three ways of pairing the jets for event. The one that pairs the 
two highest and the two lowest Jets is usually wrong since the g are mainly 
at low pt and hence is removedJlj The distribution for the two remaining 
combinations, Figure shows drop a near the endpoint for the correct com- 
bination of jets, shown as a dashed curve. The error on the mass was calculated 
by varying Mg, generating a new distribution, and using a Komogorov-Smirnov 
test to find whether the two distributions could be distinguished. The result 
of this analysis gives an estimated error 

M s - M-±/Mx° 2 = 434^ 5 16 ° ± 4.5 GeV . 

(The Komogorov-Smirnov test is a standard statistical test to see whether two 
distributions are identical. For a derivation and discussion, see KnuthEd) 

7 Global Fits to Determine SUSY Parameters 

Section |^ has described a number of precision measurements that can be made 
at the LHC by relating features of kinematic distributions to combinations 
of SUSY masses. Quite a few other such measurements have been developed 
for these points, and many more would surely be found with the incentive of 
real data indicating the discovery of SUSY. It will also be possible to measure 
many other distributions that cannot be related to combinations of masses 
in such a precise way but which can be used to constrain the SUSY model. 
Given experimental observation of such signals, one would certainly use all 
the available data in a global fit to determine the SUSY model, just as the 
parameters of the Standard Model have been determined by global fits at LEP 
and SLC. 

Making such a fit requires generating large samples of many possible SUSY 
signals, combining them with the Standard Model background, and comparing 
the results with the distributions for the assumed SUSY signal. This requires 
far more effort than can be justified at present. Instead, a much simpler ap- 
proach has been adopted. Samples of signal events have been generated for 
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Figure 36j — lErrors on SUGRA parameters from the Hinchliffc et al.0 (circles) and 
Froidevauxr 5 ! (squares) fits. 



each of the five LHC points. For each point, precision measurements of com- 
binations of masses have been studied, and the errors on these combinations 
have been estimated. Then, SUGRA models have been generated choosing 
the parameters at random, the masses have been calculated, and the distribu- 
tions of the model parameters have been determined for those models giving 
combinations of masses consistent with the precision measurements. 

Two such restricted global fits have been performed. The first_fit is based 
on the set of precision measurements studied by Hinchliffe et al.cil and uses 
statistical errors appropriate to lOfb -1 , plus an assumed 3 Ge\J_theoretical 
error on the light Higgs mass. The second fit, done by FroidevauxEj as part of 
the SUSY studies of the ATLAS Collaboration, is based on all of the Hinchliffc 
et al. measurements plus a number of additional ones and uses experimen- 
tal errors scaled to three years at the full LHC design integrated luminosity, 
300 lb -1 . Froidevaux also assumes an error on the light Higgs mass set by the 
experimental error of 0.2 GeV. That is, the Hinchliffe et al. fit is conservative 
both in regards to what is included and in regards to the integrated luminosity 
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The Froidevaux fit uses more distributions and the full LHC luminosity and 
it assumes a theoretical Higgs mass error much smaller than that currently 
available, so it requires a substantial improvement on the calculation of the 
Higgs effective potential in terms of the underlying SUGRA parameters. Nev- 
ertheless, even the results of Froidevaux might be improved by additional work 
and by studying actual data. 

The results of these fits are shown in Figure 36, which displays the frac- 
tional errors for the Hinchliffe et al. and Froidevaux fits on each of the SUGRA 
parameters except A$ at each of the five LHC points. Both fits do very well on 
m!/ 2 , the gaugino mass that sets the overall SUSY mass scale. The Hinchliffe, 
et al., fits do quite poorly on mo at Points 1 and 2 because the squark masses 
are insensitive to it and the slepton masses are too heavy to enter into any of 
the cascade decays. The Froidevaux fits do rather better because they include 
more information. The fits to tan /3 and sgn fj, uses Mh heavily and so are sen- 
sitive to the assumed theoretical errors, which dominate over the experimental 
ones. Finally, neither fit gives a significant constraint on Aq. It is possible 
to constrain the weak-scale A t , A},, and A T , but the renormalization group 
equations imply that these are only weakly related to Aq and instead flow to 
quasi- fixed points. 

Figure 36 shows that the LHC has the potential to determine the param- 
eters of the SUSY model with good precision, at least for the cases studied. It 
is important to realize, however, that the use of a simple model is important 
in obtaining these results. Most of the measurements depend mainly on the 
gluino and squarks and on their main decay products, the lighter gauginos, Xii 
%2) or xf, and the light Higgs boson h. It would be a useful exercise to vary 
the SUGRA model, e.g., by allowing splittings among the scalar masses, and 
to see how well these could be constrained. 



8 Example Point with Large tan/3 

When the five LHC points were selected, ISAJET and PYTHIA still neglected 
some mixing effects due to the yi, and y T couplings. Hence it was necessary to 
choose points with tan/3 < 10, for which these effects are small, since 
gmt grrib gm T 

Vt ~ \/2M w smf3 ' Vh ~ V2M w cos(3 ' Vr ~ V2M W cos /3 ' 

More recent versions of ISAJET have incorporated these mixing effects and so 
can be used for all values of tan/3. Large tan/3 has significant consequences 
for the phenomenology of the minimal SUGRA modela 

• The renormalization group equations produce smaller bL,R and tl,r soft 
masses. 
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Figure 37: Branching ratios for ~ * Xi e_l ~ e (dashed) and Xi T+T including T-j^T-T (solid). 

• The squark and slepton masses have larger off-diagonal terms, increasing 
the mixing and so reducing the b\ and f\ masses. 

• The combination of smaller masses and larger Yukawa couplings enhance 
b and r decays. 

• Ma and related masses are also reduced. 

Perhaps the most important consequence of these changes is that it is possi- 
ble to have M(n) < M(x§), M(xf) but M{£) > M(x%),M(xf), leading to 
dominant r decays of the gauginos. 

"Point 6" is a minimal SUGRA point with mo = TO1/2 = 200 GeV, Aq = 0, 
tan (3 — 45, /i < 0. For this choice of mo and m 1 / 2 , W, Z and h decays of the 
X2 an d xt are kinematically forbidden. Choosing tan (3 = 45 makes the fi light 
enough that r decays of these gauginos are dominant, so that the X2 ~~ * ^i? 
and — * ^i v t decays are dominant, as can be seen from Figure 

Thepdpminant decays of t's are into leptons or into one or three charged 
hadronsO 

t~ — > e~v e v T 17.8% 
tx-v t 11.3% 
p~v T 25.2% 
aj"(1260) 16.3%. 

The ai(1260) is an isovector, axial vector meson which dominantly decays into 
pir. Thus t's can be identified experimentally as hadronic jets with a mass less 
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Figure 38: Effective mass M c fj for Point 6 signal (solid) and Standard Model background 
(dashed). 

than M(t) = 1.777 GeV and with either one charged track or three charged 
tracks with net charge ±1. 

As is normally the case at the LHC, gluino and squark production domi- 
nate at this point. The same basic cuts are made as for the points discussed 
previously: 

• > 4 jets with p TA > 100 GeV, p T .2,3A > 50 GeV. 

• Missing energy f T > 100 GeV. 

• Transverse sphericity St > 0.2. 

• M e ff = fr + Pt,i + Pt,2 + Pt.3 + Pt,4 > 500 GeV. 

After these cuts, discovery is straightforward via multiple jets plus Jfix- Fig- 
ure [58] shows the distribution of the variable M e g defined previously for the 
signal and the sum of all Standard Model backgrounds. Evidently the signal 
dominates for large M e g without using any r identification cuts. 

After the above cuts, the signal is selected by requiring M e g > 500 GeV, 
two hadronic r decays, and no other isolated leptons. The r's were selected 
using information from the event generator, since the actual detector cuts 
require rather a detailed analysis]] Background from two independent decays, 

d This is in progress. It is expected that backgrounds from misidcntificd jets are not very 
important since the SUSY signal dominates before any r cuts are made. 
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Figure 39: Visible t + t (solid) and t^t^ mass (dashed) for Point 6 and Standard Model 
background (shaded). 
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Figure 40: Visible t + t~ — r^r^ mass for Point 6 signal (solid), X2 decay contribution 
(dashed), and heavy gaugino contribution (dash-dotted). 



either from other SUSY channels or from Standard Model background, were 
removed by subtracting same-sign events from opposite-sign ones. Finally, only 
r decays into three charged particles (mainly a\v T and the 37r^ T ) were used. 
Since the 3-prong branching ratio is only about 15%, this involves a substantial 
loss in statistics, but the visible 3-prong momentum is closer to the true r 
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momentum. The opposite-sign and same-sign visible mass distributions after 
these cuts for the Point 6 signal and Standard Model background are shown in 
Figure 39. The difference of the opposite-sign and same-sign signals is shown in 
Figure 4C. This figure also shows the contributions from %!] — » tit — > Xi T+T ~ 
and from events involving heavier gauginos. The former is responsible for the 
clear edge in the experimental distribution, while the latter dominates at larger 
masses. 

Point 6 and similar points at which r decays dominate will clearly be more 
difficult to study than points with decays into electrons and muons. More 
work is needed both to optimize the strategy and to understand the detector 
implications. Nevertheless, the results presented here arc not unencouraging. 
Preliminary results indicate, however, that it will be difficult to detect any 
signals for such points at the Tevatron. 

9 Higgs Bosons at LHC 

Since the main reason for introducing SUSY at the weak scale is to cancel the 
quadratically divergent corrections to the Higgs masses, observing Higgs bosons 
at the LHC is just as important as observing SUSY particles. It is also a good 
deal more difficult, partly because the couplings to light quarks and hence the 
production cross sections are small, and partly because only some rare decay 
modes can be observed over the Standard Model backgrounds. Indeed, the 
ability to .search for Higgs bosons has been the main design criterion for the 
detectors. BB 

In ±hc Minimal Supersymmetric Standard Model, there are five Higgs 
bosons ,c3 h, H, A, and H ± . These are described at tree level by just two 
parameters, generally taken to be the mass Ma of the pseudoscalar Higgs bo- 
son and tan /?, the ratio of vacuum expectation values. The renormalization 
group equations and radiative electroweak symmetry breaking relate Ma to 
the SUGRA parameters. The other masses are then given by 
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The mixing angle a between h and H is 
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Table 6: MSSM Higgs couplings relative to a Standard Model Higgs of the same mass. 
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The couplings of the MSSM Higgs bosons relative to the Standard Model ones 
can all be expressed in terms of a and /3, as shown in Table @. Note that for 
Ma ^> Mz : a ~ (3 — ^ir, so that the light Higgs h has couplings identical to 
those of a Standard Model Higgs of the same mass. 

9.1 Observing Standard Model Higgs Bosons 

The Standard Model Higgs is simpler than the MSSM Higgs sector and has 
been studied in more detail by ATLAS and CMS, so it makes sense to discuss 
it first. The dominant production is by gg fusion through a t-quark loop, 
although WW and ZZ fusion also contribute for large masses. The same 5<r 
discovery limit discussed in Section |5.2| will be used here. The .search for the 
Standard Model Higgs naturally divides into three mass regionsocll: 

M H > 2M Z 

The dominant Higgs decay modes in this mass region are WW and ZZ, 
with tt at most ~ 10% and all other modes very small. The signal for 
H — ► ZZ — > £ + £~£ + £~ is robust, with the only important background be- 
ing the ZZ continuum. It may also be possible to use the H — > ZZ — > £ + £~vv 
and H — > ZZ — > £ + £~qq modes, although the backgrounds in these channels 
are much more significant. The only problem in this region is that a very heavy 
Standard Model Higgs is strongly coupled to longitudinal WW and ZZ pairs, 
so that 

16ttV2 

Thus at very large masses, Mji > 600-800 GeV, the Higgs becomes so broad 
that it ceases to be a well defined particle and vanishes into the WW or ZZ 
continuum. Recent precision electroweak data from LEP indicateSp-Lhat the 
Standard Model Higgs is lighter than 420 GeV with 95% confidence.ca In this 
mass region, detection of a Higgs decaying into four leptons should be easy. 
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Figure 41: Standard Model H -> ZZ* signals for M H = 130, 150, and 170 GjsY, showing 
both the average signal (top) and the signal for a typical experiment (bottom) .LJ 



130 GeV < M H < 2M Z 

In this mass region, the dominant Higgs decay mode is H — > bb, but the 
background from QCD bb production is overwhelming. The best mode for 
detecting such a Higgs is H — > ZZ* — > i + l~i + t~ . While the rates are small, 
isolation cuts can be used to reject the Zbb background. Then, after the cuts 
M12 « M z and M 34 > 10 GeV to reject the Z^* continuum, the S/B is very 
good, as can be seen in Figure although the signal is quite small. It seems 
clear that one can eventually detect a Higgs in this range, although it may 
take more than one year for Mr » 170 GeV, for which the W + W~ mode 
dominates. 

For low statistics processes such as H — * ZZ*, it is important to treat the 
statistics correctly. The right question to ask is, "What is the probability for 
the background to fluctuate up to the expected signal?" For large numbers of 
events, the answer is given by using Gaussian statistics to calculate the sig- 
nificance, S/VB, in standard deviations. For low numbers of events, however, 
it is important to use Poisson statistics and then to translate the probability 
to an equivalent Gaussian significance. Note, for example, that the Poisson 
probability for a background of 1 event to fluctuate up to 6 events is 

P(5+l|l) = ^- = 5.1x IO" 4 , 
G8 



Table 7: H — > ZZ* event numbers. 



Mass (GeV) 


120 


130 


150 


170 


180 


Signal (3 x 10 4 pb _i ) 


2.6 


12.5 


34.9 


10.2 


26.5 


it 


< 0.1 


< 0.1 


< 0.1 


< 0.1 


< 0.1 


Zbb 


0.1 


0.2 


0.2 


0.2 


0.2 


ZZ*, Z~/* 


1.6 


2.8 


3.3 


3.5 


3.3 


Background (3 x 10 4 pb _1 ) 


1.7 


3.0 


3.6 


3.7 


3.6 


S/VB (3 x lC^pb- 1 ) 


2.0 


7.2 


18.3 


5.3 


14 


Signal (10 5 pb _i ) 


5.2 


24.5 


68.5 


19.9 


51.9 


Background (10 5 pb _1 ) 


4.7 


8.2 


10.0 


9.5 


9.0 


S/y/B (10 5 pb _1 ) 


2.4 


8.5 


21.7 


6.5 


17.3 



compared to 5<r Gaussian probability of 5.7 x 10~ 5 . The probably that an 
expected signal S + B = 6 fluctuates down to the background B < 1 is much 
larger: 

P(0|6) + P(l|6) = e~ 6 [l + ^] = 0.017 . 

Since there are many possible fake signals but presumably only one real one, 
this is acceptable. 

The Gaussian significances shown in Table are thus slightly too opti- 
mistic, but nevertheless they show that it is possible to discover a Higgs boson 
decaying through ZZ* into four leptons for Mjj > 130 GeV. The maximum 
luminosity is needed for Mjj ~ 170 GeV because the WW channel is open, 
reducing the ZZ* branching ratio. The acceptance used in this table may be 
somewhat optimistic, particularly for muons, but this can be compensated by 
more running time. 

M H & 130 GeV 

Below about 130 GeV the H -> ZZ* -> £+£-£+£- branching ratio becomes 
too small to be observable, and the dominant bb and rr modes are swamped 
by background. Thus in this mass range the only way to observe the Higgs is 
via the 77 mode. This mass region is particularly important for SUSY, since 
the light Higgs boson is always in this range and generally is almost identical 
to a Standard Model Higgs of the same mass. There are large backgrounds 
from the QCD 7-7 continuum which can only be overcome with excellent 7-7 
mass resolution. There are also huge potential backgrounds from the 7-jet and 
jet-jet backgrounds. Detection of H — > 77 is an extraordinary challenge, and 
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Figure 4pJri 120 GeV Standard Model Higgs signal and 77 background with the ATLAS 
detectorH 
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Figure 43: Background-subtracted signals for H 1—1 77 with the CMS detector and cross 
sections required for 5, 7, 10, and 15a significance.113 



the need to do so has driven the design of the electromagnetic calorimeter — 
and, to a lesser extent, of the whole detector — for both ATLAS and CMS. 

The first problem in detecting H — > 77 is to overcome the 77 continuum. 
The lowest order QCD process is qq — > 77, but it is important to include 
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the (gauge invariant) higher order process gg — > 77 through a box graph be- 
cause g(x) is very large at small x. The background calculation also includes 
7 bremsstrahlung from quarks in the leading-log approximation. The follow- 
ing cuts are made in the ATLAS analysis to reduce this real 77 background 
compared to the signal: 

• PT ,i > 40GeV, pt.2 > 25GeV, < 2.5. 

• Pt,i/(pt,i +Pt,2) < 0.7 to minimize bremsstrahlung backgrounds. 

• Exclude regions of the calorimeter with degraded resolution. 

• Choose AM for each mass to maximize S/^/B. 

With these cuts and with no fake 77 background, even the most favorable 
mass, Mjj = 120 GeV is hard, as can be seen from Figure |4^. Lower masses, 
Mjj < 100 GeV, are the most difficult, as can be seen from Figure [l3|. For a 
Standard Model Higgs in this mass range, one or more years at full luminosity 
are required to reach a 5a significance. For all masses the S/B ratio is poor, 
but because the peak is narrow, the side bands can be used to determine 
the background. Thus, the search limits are set by statistics and not by the 
understanding of the background. 

The potential jet-jet and 7-jet backgrounds are very large; a 7/jet rejection 
of ~ 10 4 is essential to reduce these to a fraction of the QCD 77 continuum. 
Very detailed simulations of the isolation cuts and of the rejection of single 
and multiple tt° using the detailed properties of the electromagnetic showers 
indicates that the required rejection can be achieved. An efficient trigger on 
two photons is also possible. But detecting this signal is probably the most 
difficult challenge for ATLAS and CMS. 

9.2 Observing SUSY Higgs Bosons 

In the MSSM the Higgs sector at tree level depends only on two parameters, 
Ma and tan/3. In the SUGRA model it is usually true that Ma ^S> Mz, so 
that the light Higgs h is very similar to a Standard Model Higgs of the same 
mass, while the rest of the Higgs bosons are heavy, Mjj ~ M H ± w Ma- Many 
modes might be used to searcLior the SUSY Higgs bosons. Figures |44| and 
show the 5a reach contoursEj for various signatures in the M^-tan/3 plane, 
calculated assuming that_all SUSY masses are ~ 1 TeV so that SUSY decay 
modes do not contribute.LJ Each mode has a statistical 5a significance on the 
shaded side of the corresponding curve. The plot does not show constraints 
such as the rate for b — ► S7 that limit the Higgs parameter space for the 
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Figure 44: Possible Higgs iStt discovery contours for 30 fb -1 in various modes assuming all 
SUSY masses are ~ 1 TeVc2l Each signal is observable on the shaded side of the correspond- 
ing curve. 



assumed SUSY masses. Perhaps the only thing clear from these figures is that 
many modes can be used to search for SUSY Higgs bosons. By making use of 
all of the modes, it appears to be possible to discover at least one Higgs boson 
over the entire M^-tan/3 plane by combining LEP at 195 GeV and ATLAS 
and CMS with 300 fb -1 eachJla Over much of the parameter space, more than 
one Higgs can be observed, but it is generally not possible to observe all of 
them. 

All of the Standard Model Higgs search modes will be used to search for 
SUSY Higgs. In particular, in the SUGRA model the light Higgs h is almost 
always in the range covered by the 77 mode because Ma is large, so that the 
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Figure 45: Possible Higgs fm discovery contours for 300 fb -1 in various modes assuming all 
SUSY masses are ~ 1 TeVc2l Each signal is observable on the shaded side of the correspond- 
ing curve. 



h couplings are similar to the standard model ones. Hence the search for the 
h is similar to that for a Standard Model Higgs in the same mass range. The 
5<t discovery limits are indicated in Figures and [|5]. 

For the heavy Higgs bosons H° and ^4°, the ZZ modes are suppressed 
in the MSSM compared to the Standard Model. Because of CP invariance, 
A y4 ZZ, and the H coupling to WW and ZZ is proportional to cos(/3 — a), 
which vanished for large Ma- Hence the four-lepton modes are of limited 
value for SUSY Higgs, although there are regions indicated in Figures [44| and 
|45| where they can be observed. 

Standard Model Higgs decays into rr are probably not possible to observe. 
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The H and A decays to tt are enhanced for large tan j3, as can be seen from 
Table ^, and the competing WW and ZZ modes are reduced or eliminated. 
Since the r always decays to a missing v r , it is not possible to reconstruct the 
tt mass in general. But if the Higgs is produced with high pr, pt ~ Mh, 
then the two t's are not back to back and it is possible to resolve the missing 
energy I^t along the observed r directions, reconstructing the two r momenta 
and hence the tt mass. The mass resolution for such events has been studied 
with full simulationo The dominant tt background is Z* /j* — ► tt, and this 
can be determined using Z — > ee, /i/i, and the known r decay modes. Hence, 
the statistical errors should be dominant, and the 5cr discovery limits indicated 
in Figures |4| and [l5] should be reliable. 

The process Wh — > I vbb is another one that is difficult to observe for Stan- 
dard Model Higgs, even though the lepton tag greatly reduces the Standard 
Model bb background. This signal is enhanced for large tan/3 in the MSSM 
and is at least statistically observable in the regions shown in Figures |4| and 
ff5| . But it must be noted that there is a large ti background: e.g., there are 
475 signal events on 12700 background for a 100 GeV Standard Model Higgs. 
While this gives a statistical significance S/^/B = 4.2, it is necessary to know 
the background to better than 1% in order to be sure that the signal is real. 
While the S/B for H — > 77 is also poor, the signal peak in that case is very 
narrow, so the sidebands can be used to determine the background. For bb 
decays the mass resolution will be of order 10%, and the ti background can 
produce background that varies significantly over this range. While it may be 
possible to observe this signal, the small value of S/B implies that it will be 
very difficult. 

10 Conclusions 

If SUSY exists at the electroweak scale, i.e., with squark and gluino masses less 
than 1-2 TeV, it should be straightforward to find signals for it at the LHC 
in the jets + fx channel and perhaps in many other channels. Discovery of a 
deviation from the Standard Model should be possible with an integrated lumi- 
nosity of 10fb _1 or even less for masses below 1 TeV. In many cases, it should 
be possible to determine combinations of masses from features of kinematic 
distributions, giving precision measurements of these mass combinations. If 
the SUSY model turns out to be simple, it will also be possible to determine 
its parameters from such precision measurements. Shortly after the LHC starts 
operation, either SUSY will become a central part of particle physics and of 
every subsequent TASI Summer School, or it will be relegated to an obscure 
corner of mathematical physics. 
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The LHC will mainly produce gluinos, squarks, and their main decay prod- 
ucts, the light gauginos, X11 X2; an d Xi • The dominant backgrounds for SUSY 
signatures come not from Standard Model processes but from other SUSY pro- 
cesses. For some choices of the SUSY model, it will also be possible to detect 
other SUSY particles, including some or all of the sleptons and the heavier 
gauginos. However, it is generally not possible to detect the whole SUSY and 
heavy Higgs spectrum. Thus, some of the conclusions from any LHC SUSY 
analysis will probably be model dependent. 

The Next Linear Collider — or the Next Lepton Collider if a muon collider 
should turn out to be preferable — probably will be completed much later than 
the LHC. An NLC generally can detect all SUSY particles which are kinemat- 
ically allowed. Using the fact that SUSY particles are pair produced aWhtj 
beam energy allows one to reconstruct the \i mass as well as other masses.BEil 
Such a machine may have to measure gauginos, Higgs 's, and perhaps sleptons 
with masses ~ Mg This would require y/s > 2M§. While this scale is not 
now known, it could be 2 TeV. Of course, if the SUSY model turns out to be 
complicated, a lower energy machine capable of studying the sleptons and light 
gauginos could prove essential. There are exciting experimental prospects for 
SUSY not just at the LHC but at future machines. 

I wish that I could be as optimistic about the future of particle physics in 
the United States. 
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